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Abstract. We consider the Einstein deformations of the reducible rank 
two symmetric spaces of noncompact type. If M is the product of 
any two real, complex, quaternionic or octonionic hyperbolic spaces, 
we prove that the family of nearby Einstein metrics is parametrized by 
certain new geometric structures on the Furstenberg boundary of M. 



1. Introduction 

This paper is the first in a series to investigate the deformation theory of 
Einstein metrics asymptotically modelled by Riemannian globally symmetric 
spaces of noncompact type and of arbitrary rank. In the special case of 
real hyperbolic space, and the slightly more general setting of conformally 
compact asymptotically hyperbolic manifolds, this has been the focus of 
extensive study over the last fifteen years; this attention is due both to 
the many deep connections with conformal geometry, cf. [10], as well as 
the central role that these 'Poincare-Einstein' spaces play in the AdS/CFT 
correspondence in string theory, for which the proceedings [3] provides a 
good introduction. More recently, some of this analysis has been extended to 
the other rank one noncompact symmetric spaces [2] . Some recent advances 
in linear analysis on symmetric spaces has now made it reasonable to attack 
this problem in greater generality. 

The rank one globally symmetric spaces of noncompact type are the real, 
complex, and quaternionic hyperbolic spaces, denoted MH m , CH m , MH rn , 
respectively, and the octonionic hyperbolic plane OH 2 . Each has curvature 
bounded between two negative constants and is diffeomorphic to an open 
ball B n+1 . There is a family of 'asymptotically K hyperbolic' (K = R, C, H 
or O) - or briefly, AKH - metrics, with asymptotics modelled on ¥LH m , each 
of which induces a geometric structure on the sphere at infinity, S n . In the 
real case, this geometric structure is a conformal class; in the complex and 
quaternionic case, it is a CR or quaternionic contact structure, as described 
in [2]. There is an octonionic contact structure on S 15 , but it is rigid. We 
call these boundary structures either the conformal infinity data of the AKH 
metrics, or alternately, G-conformal structures, where G is the semisimple 
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Lie group associated to ~KH m . One main result of [2] is that for each G 
conformal structure near the standard one on S n there is a unique AKH 
Einstein metric near to ~KH m , and with that conformal infinity data. This 
correspondence between AKH Einstein metrics and G-conformal structures 
can be regarded as an asymptotic boundary problem, and the assignment of 
the interior Einstein metric to the conformal infinity data a sort of nonlin- 
ear Poisson transform. There are other very interesting, and more subtle, 
problems of this type: we mention in particular LeBrun's positive frequency 
conjecture, concerning self-dual and anti-self-dual Einstein metrics in four 
dimensions, which was solved by the first author in [3]. 

It is our goal in this paper to establish a similar local deformation theory 
in the first higher rank case, namely for products of the various hyperbolic 
spaces listed above. There are a number of new and interesting features 
not encountered in the rank one case, and the details of the geometry and 
analysis are already sufficiently complicated that it has seemed reasonable 
to keep this as a separate paper. In later papers we shall treat the cases 
corresponding to more general noncompact higher rank symmetric spaces. 
This mirrors the recent developments for the linear analysis (for the scalar 
Laplacian) [3H I22J |23 125 . 

At the roughest level, the proof proceeds exactly as in the AKH setting, 
by constructing a family of approximate Einstein metrics, parametrized by a 
family of boundary structures generalizing the G-conformal structures, and 
then applying the inverse function theorem to an appropriately gauged ver- 
sion of the Einstein operator. The solution of the resulting nonlinear elliptic 
equation yields the 'near product hyperbolic' Einstein metric with the pre- 
scribed conformal infinity data. There are three main new issues in carrying 
this out for higher rank symmetric spaces. The first, purely geometric in na- 
ture, involves defining the appropriate analogue of G-conformal structures. 
This relies in turn on a choice of compactification for each of these product 
hyperbolic spaces as a manifold with corners of codimension two; the new 
boundary structures are defined on the corner. However, not every one of 
these new boundary structures can be extended to an asymptotically prod- 
uct hyperbolic metric which is also asymptotically Einstein everywhere near 
infinity. Extra hypotheses on the boundary structure must be imposed, and 
even then one must solve an extra Einstein-like equation to be able to ex- 
tend this structure from the corner to the codimension one boundary faces. 
The final issue is to attain some understanding of the mapping properties 
of the linearized gauged Einstein operator on these asymptotically prod- 
uct hyperbolic metrics so that we can apply an inverse function theorem 
argument. 

In slightly more detail, if Mj is a rank one hyperbolic space, j = 1, 2, then 
the correct notion of conformal infinity data on Mi x Mi is what we call a 
[G\ x Cr2)-conformal structure on the Furstenberg boundary dM\ x dMi- 
This is defined carefully in § [5j To construct an approximate Einstein met- 
ric associated to one of these, the first step is to extend this structure from 



A NONLINEAR POISSON TRANSFORM 



3 



the corner to the full boundary (Mi x DM2) U {dM\ x M2); this involves 
solving an auxiliary PDE on these hypersurface boundaries, which is a cou- 
pled version of the Einstein equation on each factor Mj . Once this has been 
done, we can extend the conformal infinity data to an approximate solu- 
tion of the problem, i.e. a metric which induces this designated structure on 
the Furstenberg boundary, and which is asymptotically Einstein uniformly 
near infinity. The next step is to determine the mapping properties of the 
linearized gauged Einstein operator on weighted Holder spaces in order to 
perturb this asymptotically Einstein metric to an exact one. This is done us- 
ing the geometric parametrix approach developed by the second author and 
Vasy, as referenced above, which involves an adaptation of the techniques 
of iV-body scattering theory to study elliptic theory on symmetric spaces 
of rank greater than one. A subtle but important complication is that if 
either of the factors is quaternionic, then the approximate Einstein metrics 
corresponding to different (Gi x G2)-conformal structures are not mutu- 
ally quasi- isometric, even up to diffeomorphism. Because of this, even the 
function spaces vary in a nontrivial way as we vary the boundary (G\ x G2)- 
conformal structure. This means that we must analyze the Green function 
for the linearized gauged Einstein operator not just at the exact product 
metric, but for all nearby approximate Einstein metrics. This necessitates 
that we carry out a parametrix construction at 'near product hyperbolic 
metrics', which is more complicated than doing it just at the product hy- 
perbolic space only. (This difficulty is already present in the quaternionic 
hyperbolic case [2] ; there is a way to circumvent it then which unfortunately 
does not generalize to this product setting, so the parametrix construction 
seems unavoidable here.) 

General information about the geometric analysis behind the deformation 
theory for Einstein metrics, particularly in the compact setting, can be found 
in pp. We follow a slightly different route developed in [2]. Suppose that g 
is Einstein, i.e. Ric(g) — Xg = for some real number A. If h is a sufficiently 
small symmetric 2-tensor, then to make the equation h 1— > Hic(g+h)—X(g+h) 
elliptic we supplement it with the so-called Bianchi gauge condition B 9 (h) := 
5 9 h — ^dti 9 h = 0. Equivalently, we look for solutions of the nonlinear elliptic 
system 



It is not hard to show that if A < and B 9 (h) — > at infinity, then solutions 
of this equation correspond to Einstein metrics g + h in Bianchi gauge with 
respect to g. One advantage of this gauge is that the linearization of N 9 at 
h = is the particularly simple operator 



where the final term on the right is the usual action of the full curvature 
tensor for g on symmetric 2-tensors. Throughout this paper, this operator 
will be called the linearized gauged Einstein operator. 



N 9 {h) := Ric(g + h) - X(g + h) + (5 9+h )*B 9 {h) = 0. 
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Our principal result is the 



Main Theorem. Let (Mj,gj) be an AKjH Einstein space, j = 1,2, and 
let c be the product (G\ x G2)-conformal structure on dM\ x dM2- Assume 
that is not an L 2 eigenvalue for the linearized gauged Einstein operator 
on either (Mj,gj), j = 1,2, or (Mi x M<i,g\ + g<i). Let c' be any other 
smooth (Gi x G2)-conformal structure which is sufficiently close to c in the 
<^2,o norm ^ Assume alsgo that c' satisfies the global integrability hypothesis 
in Definition A25\) . Then there is a near product hyperbolic Einstein metric 
g' with conformal infinity data c' , and moreover g' is unique amongst such 
metrics in a neighbourhood of the product metric g. 



The hypothesis on the linearized gauged Einstein operator is satisfied in 
many situations, in particular for convex cocompact quotients of IK hyper- 
bolic spaces and negatively curved AKH Einstein spaces. 

We note also that the techniques and results of this paper apply some- 
what more generally than when M is globally a product. It would not be 
too difficult to define a class of manifolds and metrics with appropriate local 
product conditions near the corners and boundary faces to which this defor- 
mation theory also applies. As a very simple example, M might be obtained 
by a compact topological perturbation from the product Mi x M2. However, 
we do not currently know any manifolds of this type which are Einstein, and 
so have not formulated our main result in this greater generality. 

One final comment about notation. We shall be using various classes 
of Holder spaces, often weighted by powers of boundary defining functions. 
We typically write c rf k ' a for Holder spaces on compact manifolds (or in any 
compact set), and Ag ,a when these spaces are defined relative to some com- 
plete metric g on an manifold. (In fact, the subscript g is replaced by some 
moniker for a general class of complete metrics with a type of prescribed 
asymptotic geometry.) 

The plan of this paper is as follows. In §[2] we review the geometry of AK.H 
hyperbolic spaces. §[3]contains a lengthy review of the geometric parametrix 
theory used to study elliptic theory on these spaces, which is called the KO- 
pseudodifferential calculus; we also establish some results here about the 
resolvent family of L 9 for such a metric. § H] reviews the analysis needed to 
carry out the deformation theory of Einstein metrics in the AKH setting. 
§[5]develops the notion of (Gi x Cr2)-conformal structures and some geomet- 
ric properties and estimates for the corresponding asymptotically product 
hyperbolic metrics. The extension of these structures to the codimension 
one boundary faces is the subject of § El The parametrix construction in 
the near product hyperbolic case is the topic of § El and finally, the brief § M 
finishes the proof of the main theorem. 
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2. Asymptotically K hyperbolic spaces 

This section reviews the geometry of rank one symmetric spaces of non- 
compact type, and of the more general class of Riemannian manifolds which 
are asymptotically modelled on these. 

Hyperbolic spaces and their conformal infinities. The noncompact 
symmetric spaces of rank one are commonly called hyperbolic spaces and 
written as ¥LH m , where K = R, C, H (the quaternions) or O (the octonions). 
Note that OH m exists only when m = 1, 2, and in fact OH 1 = RH 8 , so the 
only new space in this last family is the 16-dimensional octonionic hyperbolic 
plane OH 2 . As a homogenous space, ~KH m = G/K, where G is a real 
semisimple Lie group and K a maximal compact subgroup; more specifically, 

RH m = SOi jm /SO TO , CH m = SUi jm /U m , 

MH m = Sp 1)m /S Pl S Pm , OH 2 = Fj 20 /Spin 9 . 

These are the noncompact duals of the corresponding projective spaces 
KP m . 

Throughout this paper we write 

(1) d = dimit K, 
so that 

diniRKH™ = md:=n + l, 

where this last equality defines n in terms of m and d. 

The polar coordinate expression for the metric on WH m is 

(2) g = dr 2 +sinh 2 (r)7, 

where 7 is the standard metric on S n . For the analogous expression on 
the other hyperbolic spaces, denote by r] G S7 1 (S' ri ) (g) Im(K) the connection 
1-form of the Hopf bundle 

gd-i — „ S n 
I 

K pm-1 

and let 7 be the pullback of the standard metric on KP m_1 , regarded as a 
metric on the distribution @ = kerr/. The metric on M,H m , normalized to 
have sectional curvatures in [—4, —1], is given by 

(3) g = dr 2 + sinh 2 (r ) 7 + sinh 2 (2r ) rj 2 . 

The metric 7 on the distribution can be obtained as the limit as r -> 00 
of the family of metrics 7 r = 4e _2r g\rs n on S n ; note that this limit is finite 
only on and becomes infinite on any complementary direction. There 
is no natural 'origin', so e~ 2r and 7 are defined only up to a multiplicative 
factor. Indeed, once we are in the fully geometric setting below, it is most 
natural to take r as the distance from a large convex hypersurface, and then 
we see that 7 is only determined up to an arbitrary smooth positive factor; 
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thus only the conformal class [7] of this metric on Ql is well-defined. We call 
this asymptotic data [7]) the conformal infinity of g. 

Associated to the distribution 9) on S n is the bundle S> © (TS n /&) over 
S n . The Lie bracket on sections of Q) equals — drj, and thus induces the 
structure of a nilpotent Lie algebra on each fibre of this extended bundle 
which is isomorphic to the IK-Heisenberg algebra IK — Heis m_1 = IK m_1 © 
Im(IK). The metric 7 is compatible with dr] 6 Q% © Im(K) in the sense that 
the pair (dr], 7) defines a K-structure on @, i.e. a collection of d — 1 almost 
complex structures which are orthogonal with respect to 7, and which satisfy 
the algebraic relations of the basis elements in Im(K). 

G-conformal structures. The hyperbolic metrics (|2|) and ([3]) are the mod- 
els for more general asymptotically hyperbolic metrics (of type M, C, H or 
O). Before defining these, however, we first introduce terminology for the 
conformal infinity structures, which will also be used later in the product 
case. 

Definition 1. Fix the IK hyperbolic space ¥LH m = G/K . A G-conformal 
structure on an arbitrary manifold Y n , n = md — 1, is a codimension d — 1 
distribution £) C TY, with a conformal structure [7] on the fibres of *3 , such 
that the induced nilpotent Lie algebra 3i ®TY / 2> is isomorphic at each point 
to the "K-Heisenberg algebra, and any metric 7 G [7] is compatible with the 
IK structure on 3> . We say that the distribution *2) is of ^-contact type. 

This definition unifies several cases: 

• when IK = M, 'Ss is the entire tangent space and [7] is a conformal 
structure in the usual sense; 

• when IK = C, S> is a contact distribution in the ordinary sense; if rj is 
a contact 1-form which defines then the compatibility of a metric 
7 on ^ with dr] means that on £F one has dr](-, •) = 7(J-, •) for some 
almost complex structure J which is orthogonal with respect to any 
7 6 [7]; this is simply an almost CR structure on S>; 

• when IK = H, 2> is a "quaternionic contact structure" as defined and 
studied in [2j; it turns out that the conformal class [7] is completely 
determined by drj; 

• finally, the octonionic case is rigid; & is automatically locally iso- 
morphic to the standard distribution on the sphere 5 15 , and [7] is 
determined completely by dr]. 

It is important to note here that in the quaternionic case, even though 
each of the tangent nilpotent Lie algebras is isomorphic to the standard 
quaternion Heisenberg algebra, the distribution £F is not locally diffeomor- 
phic to the model structure on the sphere (unless it is standard everywhere). 
In other words, there is no direct analogue of Darboux's theorem for quater- 
nionic contact structures, and the infinitesimal equivalence of these struc- 
tures at each point does not imply their local equivalence. 
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AKH metrics. Let M be a manifold with boundary Y = dM admitting a 
G-contact structure (associated to ~KH m = G/K). We give two equivalent 
definitions of the class of complete metrics on the interior of M which induce 
a G-conformal structure on Y . 

The first mimics the polar coordinate definition of the model case. 

Definition 2. A metric g on M n+1 is called asymptotically K hyperbolic 
(or AKH for short) if the following conditions are satisfied: there is a neigh- 
bourhood % ofYin M, a diffeomorphism identifying with (1, oo) r x Y , 
and a G-conformal structure (St, [7]) on Y , such that, fixing a representative 
(jlil) of the G-conformal structure and defining 

9o(l, v) = dr 2 + sinh 2 (r) 7 + sinh 2 (2r) rj 2 

in , we have 

9 = go(i,v) + k , 

where k is in the weighted geometric Holder space e~ ur A 2 ' a for some v > 0. 
(The derivatives and norms are taken with respect to 5 f o(7i ? ?)-j 

This definition does not depend on the choice of (77, 7) in the confor- 
mal class since, replacing (77,7) by (fr],f~f) for some / £ ^°°(M), / > 0, 
changes the model (up to diffeomorphism) by an error which is G{e~ r \ The 
pair ([7], 7/) (or more properly, the triple [7], 77)), is called the conformal 
infinity of g. 

The alternate definition simply replaces the radial variable r by x = e~ r , 
which is a defining function for Y in M (recall, this means that x > in M, 
x = only on Y and dx 7^ there); thus 

glx 2 + 7 T] 2 

(4) 9o(l, V) = — — 2 — + ^4' 9 = 9o(l, V) +k, 

where k £ x u h? ,a . (Again, norms and derivatives are with respect to 
50(7, T]).) This will be more useful from our point of view since the boundary 
Y appears explicitly as the hypersurface {x = 0}. Near the boundary, the 
volume form of g has the form 

(5) dV g = x~ n - d dxdA Y {x) 

where dAy(x) is a family of volume forms on Y depending smoothly on x. 

A straightforward calculation, cf. [2], shows that an AICH metric has 
curvature tensor which is asymptotic to that of ~KH m to order ff(e~ ur ) = 
6(x v ). 

The complex hyperbolic metric on the unit ball B in C m has a slightly 
different form in standard Euclidean coordinates, and it is worth explaining 
the difference. This metric has Kahler form 

ddp dp dp 
P P P ' 

where p = ^(1 — \z\ 2 ) is a defining function for dB. As a Hermitian metric, 
the first term blows up only like 1/p; its leading coefficient is the Levi form, 
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which is positive definite on Q. The second term, which blows up at the 
faster rate l/p 2 , vanishes on (the radial extension of) @, and is positive on 
the directions spanned by dp and dp, or equivalently, on the span of d p and 
id p . There are analogous expressions for the Bergman and Kahler-Einstein 
metrics on any strictly pseudoconvex domain. 

The obvious discrepancy with (j3|) is resolved by setting x = yfp. This 
accords with the fact that the geodesic distance function r for the hyperbolic 
metric is comparable to — \\ogp rather than — log p. More bluntly, the 
standard structure on the closure of the Euclidean ball (or any strictly 
pseudoconvex domain) induced from its inclusion in C m is not quite the 
right one for our purposes. 

3. Linear elliptic theory on asymptotically hyperbolic spaces 

We now describe the structure of the Green function and mapping prop- 
erties for the linearized gauged Einstein operator on an asymptotically K 
hyperbolic space. 

There are several ways to approach linear elliptic problems of this type. 
Because the underlying geometric structure is asymptotically rank one, cer- 
tain features of the operators in question are dominated by their radial be- 
haviour, which is one-dimensional and hence more readily tractable. Using 
this, the first author [2] carried out a detailed ODE analysis for the radial 
part of the relevant operators on each K hyperbolic space to capture the de- 
cay of the corresponding Green functions, from which the required mapping 
properties can be deduced. In principle, the same general ideas should work 
for higher rank geometries, but the radial parts of these operators are then 
multi-dimensional and must be studied differently This will be done later 
in this paper via the more general techniques of geometric microlocal analy- 
sis. We review these methods in the setting of asymptotically IK hyperbolic 
geometry, even though simpler methods are available there, because this is 
a good warm-up for the construction in the product case below, but also 
since we require certain more subtle estimates on the resolvent family here 
which are used in the product analysis of §7. 

Before embarking on all of this, let us say a few words about the general 
strategy. Local elliptic theory, or elliptic theory on compact manifolds, can 
be developed entirely via Schauder estimates. A satisfactory understanding 
of the global mapping properties for a Laplace-type operator L on a com- 
plete noncompact manifold requires not only this local theory but also some 
information about 'far-field' effects. Roughly speaking, one needs estimates 
at infinity for solutions of Lu = f, even when / £ ^q°- Supposing for 
simplicity that L is actually invertible on L 2 , then its inverse is represented 
by an integral operator / u(z) = J G(z, z')f(z') dz' . This integral ker- 
nel, G, is a distribution on M x M, and is called the Green function for 
L. Its structure near the diagonal {z = z'} is exactly the same as in the 
compact case, but the interesting mapping properties of L are determined 



A NONLINEAR POISSON TRANSFORM 



9 



by its asymptotics as (z, z') — ► oo in any direction in M x M. For example, 
for nonlinear problems one must usually understand the invertibility of L on 
Holder spaces rather than Sobolev spaces. It is not so easy to deduce Holder 
boundedness from I? boundedness on a noncompact space, but this follows 
from pointwise estimates for the off-diagonal asymptotics of G. Thus it is a 
fundamental goal to determine these asymptotics, at least for special classes 
of complete spaces. 

The existence of G may be known by abstract or indirect methods, e.g. 
using Hilbert space theory and a Bochner-type argument, but these usu- 
ally give little information about the asymptotic structure of this Schwartz 
kernel. A parametrix for L is an approximation to G, or slightly more gener- 
ally (in case L is only Fredholm) an operator which inverts L up to compact 
errors (though note that compactness depends on the function spaces on 
which these act). A parametrix construction produces operators for which 
we have good pointwise control on the Schwartz kernels. In these special 
geometric settings, such a construction proceeds by solving a sequence of 
model problems, and using their solutions to construct successively better 
approximations to the true (putative) Green function. To organize this 
information in a useful manner, we begin by defining a certain compacti- 
fication M 2 of M x M. This compactification is a manifold with corners, 
and the model problems appear as the induced operators on its different 
boundary hypersurfaces. The fine pointwise structure of a parametrix is 
encoded in the statement that it is a conormal (or even better, a polyhomo- 
geneous conormal) distribution on M 2 . A posteriori one also deduces this 
same regularity structure for the Green function itself. The technicalities of 
this construction involve defining a pseudodifferential calculus on M which 
is large enough so that one may carry out some sort of parametrix construc- 
tion for 'fully elliptic' operators. These operators are characterized by the 
regularity properties of their Schwartz kernels on M 2 , and the main work 
consists in verifying the usual properties, i.e., composition, boundedness, 
etc. This can be done in a number of situations under fairly strong assump- 
tions on the ambient geometry, all satisfied in this asymptotically hyperbolic 
(or higher rank symmetric) setting. In this section we explain how to carry 
this out for the asymptotically K hyperbolic geometries. 

We assume in this section that L is a generalized Laplacian acting between 
sections of bundles E and F over M which are associated to the bundle of 
orthonormal frames (or else as a parallel subbundle when the holonomy of 
the metric is reduced) via a representation of the orthogonal group; such 
bundles come equipped with the Levi-Civita connection V. For example, 
these could be parallel subbundles of the full tensor bundle. A geometric 
differential operator is a linear combination of powers V fc of the covariant 
derivative on one of these bundles, with coefficients determined by the metric 
and curvature tensor. In particular a generalized Laplacian on a geometric 
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bundle E is an operator of the form 

L = V*V + R, 

where R is a symmetric endomorphism on E constructed from the curvature 
tensor associated to V. (However, all that we say here adapts easily to first 
order Dirac-type operators, and to many other operators besides.) In this 
section we describe the construction of a parametrix for L, and explain how 
it gives information about the Green function and mapping properties on 
various function spaces. At the end of this section we also collect some 
additional facts about the Schwartz kernel of the resolvent family R(X) = 
(L — A) ; this is an important ingredient in the construction of the Green 
function in the product case. 

We conclude with some historical remarks. The analysis of elliptic uni- 
formly degenerate operators (which is the case K = M) appeared in [18] and 
|19j . but see the monograph [T7] for a more comprehensive treatment and 
many applications. The complex case was developed initially by Epstein, 
Melrose and Mendoza [9 J ; further development and ramifications of this the- 
ory are contained in the unpublished manuscript [8j. The quaternionic and 
octonionic cases have not been written down explicitly before, though the 
way to do so had certainly been clear from [9]. These various pseudodiffer- 
ential calculi are quite similar to one another, with only minor and obvious 
modifications needed between them. One small issue, which presents only 
minor difficulties, is the lack of a Darboux theorem when K = H. Using 
language introduced below, the point is that these constructions depend on 
the infinitesimal, rather than the local, identifications of manifolds with KG 
structure and the K-Heisenberg models. 

KG structures. We first describe the notion of a KG structure on a man- 
ifold with boundary M, and its ancillaries: the KG tangent, cotangent and 
tensor bundles, and the classes of KG metrics and KG differential opera- 
tors, the latter of which contains all geometric elliptic operators for any KG 
metric as elliptic elements. 

Suppose that Y = dM carries a distribution Q> of type K; let ([7],??) be 
any associated conformal infinity. Choosing an identification of a neighbour- 
hood % of Y in M as a product Y x [0, l) x , we extend this data to Hf, and 
hence write down the model AKH metric qq = go (7 , rj) as in (jl]). Now define 
the space of all smooth vector fields V on the (closed) manifold M such 
that go(V, V) is smooth (in fy) up to x = 0. This is independent of 7 and 
77, but depends on Q (and the 1-jet of its extension to the interior). 

It is helpful to write this out in a local frame. First choose a local frame 
Yy, . . . , Yi, £ = m(d — 1), for Q and another set of independent vector fields 
Zi, . . . , Z r , r = d — 1, which are complementary to S> at each point and 
tangent to each Y x {x}; the vector field d x completes this to a full basis of 
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sections of TM. Then V E ii^o if and only if 

l r 

(6) V = a xd x + bj xYj + x 2 Z k , 

j=i k=i 

where a, bj, c k are all up to the boundary. Hence is the span over 
°°(Af) of {xd x ,xYx, xY e , x 2 Z l , x 2 Z r }. 

The terminology '0-structure' comes from [9], where denotes a nonva- 
nishing section of Tq M M <g> Im(K), the pullback of which to T*dM equals rj. 
Hereafter, we let denote not only this form on dM, but also some choice 
of smooth extension to the interior. All of the notions here can be defined in 
terms of this form, or equivalently, the corresponding oriented (d — l)-plane 
bundle S of 1-forms on M at dM. Thus, for example, elements of are 
also characterized as vector fields which are smooth on M and which satisfy 
Q(V) = &(x 2 ). 

Note that Ykq is closed under Lie bracket. Next, there is a vector bundle, 
K& TM, for which ~f%Q is the entire space of smooth sections. The fibres are 
defined by 

Ke T p M = W^kg, 

where J? p denotes the space of smooth functions on M vanishing at p. (For 
a more prosaic definition, we take the sections xd x , xY{,x 2 Zj as a local basis 
of sections of Ke TM.) This bundle is naturally isomorphic to TM over the 
interior, but the natural bundle map 

K@ TM — ► TM, 

defined via evaluation, V >— > V(p), is the zero map when p G Y. The subbun- 
dle over dM spanned by {xYi,x 2 Zj} can also shown to have an invariant 
definition, and we denote (with a slight abuse of notation) by Ke T5M. 
Third, when p £ Y, the subspace J ! p "V%b is an ideal in with respect 
to bracket of vector fields; hence for such p, K& T p M is a Lie algebra, where 
the Lie bracket of two elements given as the equivalence class of the vector 
field bracket of representatives of the two individual classes, and as such 
is isomorphic to the solvable homogeneous extension Ky m of KHeis m_1 , 
the K-Heisenberg algebra. Note that K& T p dM is a nilpotent subalgebra, 
isomorphic to KHeis m_1 itself. 

In the simplest case, when K = R, this space of vector fields is usually 
called the space of uniformly degenerate vector fields, denoted and con- 
sists of all smooth vector fields on M vanishing at Y. There is no form 
now, so to avoid complicating the presentation we shall mostly discuss only 
the other cases, save for a few passing comments about the real case. The 
other familiar case is when M is a strictly pseudoconvex domain in C m ; 
the C0 structure is the CR structure on the boundary and the canonical 
Bergman or Kahler-Einstein metrics are C0 metrics (admittedly only poly- 
homogeneous rather than c € oa ). As explained earlier, one needs to take x 
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as the square root of the Euclidean distance to the boundary in order to fit 
this into the present framework. 

The dual of the K0 tangent bundle is denoted K ®T*M. Note that smooth 
sections of this KG cotangent bundle are singular in the ordinary sense: in 
terms of the dual basis of one- forms dx, Y* and Z*, 

where a,bi,Cj £ ^°°(M). Similar remarks apply to all other tensor bundles 
too. Note in particular that an AKH metric g is a section of S 2 ( K ®T*M) 
which is positive definite on K ®TM. If E is any bundle constructed func- 
torially from TM, then applying the same functorial operations to K@ TM 
yields a bundle which we denote Ke E. 

Definition 3. Let M be a compact manifold with boundary and £F a dis- 
tribution of type K on Y = dM. The space DifT^@(M) 0/K6 differential 
operators on M consists of all operators which can be locally expressed as a 
finite sum of products of elements of Ykq- If E,F are any vector bundles 
over M , then a ICO operator acting between sections of E and F is one 
which has this form with respect to any local trivialization. 

It is important for ^0 to be closed under Lie bracket for this space of 
operators to be well-defined. 

In our applications, the bundles E and F are geometric bundles, and the 
operator is a K0 differential operator between K@ E and K6 F. 

Theorem 4. Let g be any AKH metric on M , and L a geometric elliptic 
operator of order fj, between sections of two geometric bundles E and F. 
Then L £ Diff£ e (M; Ke E, Ke F). 

This result is tautological once one checks that the Levi-Civita connection 
V satisfies 

V : <*f°°(M; E) — ► ^°°(M; E ® K@ T*M). 
We leave details to the reader. 

There is a principal symbol mapping for I^k® operators, defined formally 
replacing xd x , xYi and x 2 Zj, respectively, by linear coordinates i,t]i,C,j- 
Thus, 

Diff^ (M;£,F) 3P= £ a ]a p{z){xd x Y (xY) a {x 2 Zf 

i+H+|/3|<M 

^ Ke a,(P)(z^,r,,0 = ^(^VC*. 

j+\a\+\f3\=H 

The usual calculation shows that this is a well-defined smooth function on 
K0 T*M (with values in Hom(.E, F)), homogeneous of degree \i on the fibres. 

Definition 5. The operator P G Diff^ Q (M; E, F) is called (KQ) elliptic if 
Ke a fl (P)(z; £, 77, is an invertible endomorphism whenever (£,r],Q 7^ 0. 
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Parabolic dilations and model operators. The key to the analysis of 
KQ operators is their approximate dilation invariance. More precisely, for 
any p G Y, one may define an equivalence class of dilations based at p. When 
K = R, these are ordinary radial dilations, but in the other cases the dilations 
are 'parabolic'. Using these we can define for any P G Diff^e the normal 
operator N p (P); this is a finite dimensional reduction in that it is a left- 
invariant operator on the solvable group K^ m , depending parametrically 
on p G Y. Its invertibility (for all p) is the other key hypothesis, besides 
symbol ellipticity, needed to prove that P is Fredholm. 

We begin by defining these families of dilations. The situation is simplest 
when IK = R; in this case, choose a diffeomorphism of a neighbourhood of p 
in M with a half-ball around the origin in the half-space R™ = {(s, u) : s > 
0,u£ R™- 1 }. Now use this identification and the ordinary dilation operator 
Ms : (s,u) i ^ (Ss, 5u) to define the sequence of pushforwards of the vector 
field V £%: 

(7) limfM,- 1 )^ = N p (V). 

o— >U 

From the local coordinate description of V, we see readily that N p (V) is 
a left-invariant operator defined on R™ = RJ^ n . More generally, for any 
uniformly degenerate differential operator P, define 

(8) limtM^P = N p (P); 

o— >U 

This is in the universal enveloping algebra of R^ n , and is well-defined up 
to the action of an element A G R^ n . 

In the other cases we begin by recalling the parabolic dilation structure 
on K r y m . To define this, recall that KJ? = R+ x KHeis" 1 ^ 1 (this is just the 
A- N part of the G = KAN decomposition). Choose a system of coordinates 
(s, a, u) where s > 0, a G R d_1 = Im K and u G R^™" 1 ) = K m ~ l so that (in 
K coordinates) 

(9) 6 = da + ^lm(du ■ u) 

defines the standard IK contact structure on KHeis m_1 . Thus, 

(10) sd s , s(d Ul - \lm(d Ui u)) , s 2 d aj 

is a basis of left-invariant vector fields, where in Im(d Ui u) we identify the 
vector d Ui with a vector with IK coordinates in lK m_1 , and the imaginary 
quaternions with vertical vectors d aj . The dilation is then given by 

Ms(s,a,u) = (5s,5 2 a,5u). 

(Note that these vector fields are homogeneous of degree with respect to 
Mg.) When IK = C or O, we can choose a diffeomorphism as before which 
identifies a neighbourhood of p in M with a neighbourhood of in IK5^, 
which carries the distribution @ to the model distribution on KHeis" 1-1 ; 
the model 0q above is a suitable choice for 6 on M. In the complex case 
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this uses the Darboux theorem, while in the octonion case this follows from 
the local rigidity of octonion contact structures (so that (Y, 3)) is locally 
identified with the model geometry). In terms of this identification, we 
define N p (P) by the same formula as above, arriving at an operator which 
is left-invariant on K=y and well-defined up to translation by an element of 
this group. In the last case, K = H, one has no longer the Darboux theorem 
or rigidity, but the following Lemma is proved in the Appendix. 

Lemma 6. For any quaternionic contact structure on y 4m_1 ) a nd any point 
p G Y, there exist local coordinates (o~,u) G Im(H) x H m_1 , such that the 
quaternionic distribution is given by the kernel of a 1-form © with values in 
Im(H), and the difference with the standard form 0o of the Heisenberg group 
near the origin (a = 0, u = 0) satisfies the estimate |0 — ©o| = 0(|u| 2 + |<r|). 
The two nilpotent algebra structures coincide at the point p. 

This result constructs a diffeomorphism from a neighbourhood of p to 
a neighbourhood of in M,y so that the distributions agree at the origin. 
Letting {1^} and {Y®} be local frames for *3 and the model distribution (in 
HHeis" 1-1 , extended to MJ^" 1 and then transfered to this neighbourhood), 
then clearly Yi = Y® -\-0(\u\ 2 + |<t|). It follows that the limit of the parabolic 
dilations of P is still a left-invariant operator on M,y. 

To express this more concretely, fix a boundary defining function and 
smooth vector fields Y{, i = 1, . . . , d(m — 1), Zj, j = 1, . . . , d — 1, such that 
the Yi span the extension of 3 and the Zj span a subspace complementary 
to 2> at each point. Using the obvious multi-index notation, write 

P= E a jaP {w){xd x y{xY) a {x 2 Zf, 

j+\a\+\fi\<m 

where the coefficients are assumed to be (up to the boundary). The 
values of xd x , xYi and x 2 Zj at p fix an isomorphism of K& T*M and K^, 
and 

N P (P) = a ja/3 (p)(sd s y( S Y ns 2 Z f 
j+\a>\+\P\<m 

where s G M + , and {Y®, . . . , Y^ m _ 1 ^ ) , Zf, . . . , Z®_ ± } are a fixed basis of left- 
invariant vector fields on KHeis m_1 . 

The following result is well-known in the real and complex cases. It is 
obvious in the octonionic case, and a direct consequence of Lemma [6] in the 
quaternionic case. 

Proposition 7. Let g be an AKH metric on M , and L be a generalized 
Laplace operator on M. Then, at each point of DM, the normal operator of 
L is the corresponding operator on the hyperbolic space ~KH m . In particular, 
it does not depend on the point of DM. 
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For P as above, there is a simpler family of model ordinary differential 
operators on R + called the indicial family, defined by the expression 

I P (P) = ]T a j00 (p)(sd s y . 

j<rn 

The coefficients are endomorphisms of E p . Since this is a constant coefficient 
Fuchsian operator, it is equivalent by Mellin transform to multiplication by 
a (matrix-valued) polynomial 

I P (P;C) = E a ^o(p)(O j . 

j<m 

A number ( G C, is called an indicial root if I p (P;() is singular. This is 
equivalent to the requirement that 

P(x^v(y)) = &{x^ +1 ) Vi)er(F). 

These indicial roots are fundamental invariants of P. 

Blowups and the ICG double space. There is a more sophisticated way 
to interpret the parabolic dilations, leading to a more obviously invariant 
definition of normal operators. The idea is to introduce a resolution, or 
blowup, of the product space M x M which reflects the scaling invariance 
properties of ICQ differential operators near the boundary. This provides the 
means to define the K0 pseudodifferential operators. We describe this now. 

As usual, we start with the simplest case K = R. The distributional 
Schwartz kernels of pseudodifferential operators are singular along the diag- 
onal in M x M. Unfortunately, this diagonal intersects the corner dM x dM, 
making it difficult to describe the precise structure of its singularity near this 
intersection. To remedy this we introduce a new space 

Ml = [M x M;diag aMx9M ], 

where this notation on the right indicates that we blow up M x M at the 
boundary of the diagonal. This amounts to replacing this submanifold by the 
space of inward pointing unit normal vectors; the space Mq is endowed with 
the smallest 'rf 00 structure containing the lifts of all smooth functions on M x 
M and polar coordinates around diag 9M . Thus Mq has three hypersurface 
boundaries, B\q and Bqi, the left and right faces, which are the ones lifted 
from the two hypersurface boundaries in M x M, and the new front face 
Bn created in this blowup, which is often also denoted ff. The blowdown 
map (3 : Mq — ► M 2 is a smooth mapping of manifold with corners. 

The front face Bn fibres over di&g@ M , with fibre at p G Y the set of unit 
inner normal vectors at that point; this is a quarter-sphere, the interior of 
which carries a natural projective structure. Let (x,y) and (x',y') denote 
coordinates on the two copies of M in M 2 ; we are blowing up the submanifold 
x = x' = 0, y = y', and so it is legitimate to introduce the new singular 
coordinate system s = x/x', u = (y — y')/x', x' , y' . The coordinates (s,u) 
are then projective coordinates on this quarter sphere. The normal operator 
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of P is the restriction to the fibres of B\\ of the lift of P from the left factor 
of M to M 2 and then to Mq . Thus, as in the previous definition, each N P (P) 
acts on a half-space Wf = ]R + tx R n_1 . The underlying dilation structure is 
implicit here since we are taking the normal blowup, which involves ordinary 
homothetic scaling in the tangent spaces. 

There is a similar development for the other cases, but the normal blowup 
of the boundary of the diagonal must be replaced by a blowup of this sub- 
manifold which respects the underlying parabolic dilation structure. Now, 
instead of ordinary spherical normal vectors, we use equivalence classes of 
paths converging to p, where the equivalence relationship is governed by the 
form KB. 

First define ly to consist of the smooth functions on M vanishing on Y 
and, recalling the bundle S C Tg M M determined by 0, let Is denote the 
subset of those elements / G ly such that df\ Y are sections of S. Next, fix 
p G Y and define the set of 5-parabolic curves at p to consist of those smooth 
functions 7 : [0, 1] -► M, 7(0) = p and /( 7 (t)) = for all / G I s . We 

define an equivalence relation on such curves: 

71 ~ 72 ^ /(7iW) " /(72W) = 0(t 2 ) V/ G ly; 

/(7iW)-/(72(i)) = ^ 3 ) V/€ J S . 

The space of equivalence classes is the set of inward-pointing S'-parabolic nor- 
mal vectors to Y at p, which we denote y p , and these fit together to form a 
bundle 5? over Y. Each 5^ p has a natural M + and additive structure, defined 
by 5 [ 7 (t)] = [ 7 (*t)] and [7] = [ 7l ] + [ 72 ] if /( 7 (t)) - (/( 7 i(t)) + /(72W) = 
ff(t 2 ) for all / G /y, and <^(i 3 ) for all / G Is, respectively. These do not 
define a linear structure, however, since the scalar action does not distrib- 
ute over addition. However, directly from a local coordinate calculation one 
finds that = KJ^, and scalar multiplication corresponds to parabolic 
dilation. 

Conormal distributions. We make a small diversion from the main thread 
of this section to define various classes of conormal distributions on a man- 
ifold with corners; these will appear in many places below. Let X be a 
manifold with corners and {-Hj"}vLi an enumeration of the boundary hyper- 
surfaces of X. We assume that each Hj is a smooth embedded submanifold 
with corners in X, so we can fix a global defining function pj for that face, 
i.e. pj > 0, Hj = {pj = 0} and dpj / there. 

The space i / b(X) of b- vector fields on X consists of all smooth vector 
fields V which are arbitrary in the interior of X and which lie tangent to 
all boundary faces, and hence all corners. It is helpful to describe these 
in local coordinates. Any point p G dX lies in a corner of codimension k, 
and hence is in the intersection of k faces Hj 1 , . . . , Hj k . There are local 
coordinates (x±, . . . , Xk,y) based at p with Xi = p^, i = 1, . . . , k, and where 
y = (yi, . . . , y n -k) lies in an open neighbourhood of in IR™^. In terms of 
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these, any V G % can be written as 

V = y^ j aij(x,y)x i d Xj + y^ j b i (x,y)d ye , aij,b E G < ^°°. 

In other words, 'tb(X) is spanned over by the local sections Xid Xj , 

dy e , i, j = 1, . . . , k, I = 1, . . . , n k. 

Now define the space of conormal functions of order on X 

s/°(X) = {u : Vi, . . . , V r u G L°° (X) V r > and Vj G %}. 

Any u G .e/ is smooth in the interior of X, and is in some sense tangentially 
regular at all boundaries (though it need not have a well-defined boundary 
value). Defining regularity using these 6-derivatives rather than ordinary 
derivatives allows functions like x 7 , Re(7) > 0, and |logx| s , s £ C. Next, 
for any A;-tuple a = (oi, . . . , ct^) 6 R k , set p a = p^ 1 . . . p a k k , and define 

#/ a (X) = {u = p a v : v G 

This is space of conormal functions of (multi)weight a. We also write 
st{X) =U a ^ a (X). 

A particularly useful subclass of conormal functions is comprised by the 
poly homogeneous functions. By definition, u £ si (X) is polyhomogeneous 
if near any point p on a corner of codimension k in X, u has an asymptotic 
expansion 

u ~ ^a 7/ (y)x 7 (logx/, a 7/ G < *ff°°. 

Here 7 = ( 7 i,..., 7fc ) G C fc , £ = (£ 1: . . . , 4), 

x 7 = x^ 1 . . .xf, (logx/ = | logxil^ 1 ... | logx fc | 4 , 

and where (7, i) vary over a discrete set of values in C fc x N fe which has finite 
intersection with each sector U^ =1 {7 : Re(7j) < Cj} x N fc . 

We can also talk about conormality or polyhomogeneity of a function at 
an interior p-submanifold Y C X, which by definition is locally of the form 
{x^ = . . . = Xj. = 0, y ri = . . . = y Ts = 0} (the p means 'product', i.e. Y is 
locally of product form in X). We reduce to the situation above by blowing 
up X along Y, so that u is conormal or polyhomogeneous at Y if it has 
either of these properties at the boundary face of [X; Y] which covers Y. 

Although most of the specific functions we deal with below are actually 
polyhomogeneous, we shall not emphasize or need this extra structure, ex- 
cept in the following weak sense. Fix a boundary face H and suppose that 
the weight &h corresponding to this face is 0. Let us say that u G s/g(X) 
if u G £/ a (X) as before, and that near H, u = uq + v where uq,pJj £ v G s/ a 
and uo is smooth up to H. This is simply a way of saying that u decomposes 
into a 'leading coefficient' no, which is a smooth function on H (conormal 
at all boundaries of H) and a conormal remainder term v which vanishes to 
some positive order. If ,Jtf? is a subset of the set of all boundary faces, then 
£/^(X) consists of functions with this type of decomposition at each face 
H G Jf. 
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Finally, let Y be an interior p-submanifold and Jt? some subset of the 
boundary faces, all elements of which intersect Y at dX. As before, suppose 
that on = for all H £ Jif. Define £/£,(X;Y) to consist of the functions 
u which can be decomposed as a sum u' + u" where u' £ £/^(X) and u" is 
supported in a small neighbourhood of Y and polyhomogeneous on Y, and 
such that this polyhomogeneous singularity of u" is smoothly extendible 
across all boundaries of Y. 

All of these definitions generalize immediately if u is a section of some 
smooth vector bundle over X. 

KQ pseudodifferential operators. The K0 double space provides the 
geometric setting for Schwartz kernels of KQ pseudodifferential operators. 

Definition 8. For any fx G M and set of weights a = (o"io, 001, 0) correspond- 
ing to the boundary faces Bm, An and Bn of M^q, the space ^q(M) of 
conormal KQ pseudodifferential operators on M consists of all those opera- 
tors A on M with the following properties: 

• the Schwartz kernel Ka of A is the pushforward ( under the blowdown 
(3 : Mj| e — > M 2 ) of a distribution ka on M^ Q ; 

• ka is a distribution on M^q which is conormal with respect to all 
boundaries, partially homogeneous with respect to the front face, and 
which has a polyhomogeneous singularity of pseudodifferential order 
H along the lifted diagonal diag K @, i.e. 

K A € ^ n (M]| e ,diag K e). 

Slightly more generally, we also define ^q(M) in an analogous way when 
°~ = (ciO) °~oii o"n) an d o~\\ > (however, dropping the partial polyhomo- 
geneity at Bn and only requiring that ka is conormal and vanishes to order 
an at that face). 

The action of Ka on a function / on M requires the choice of a density 7 
on M against which to integrate, so that (Ka/)(z) = J M Ka(z, z')f(z') j{z'). 
It is purely a matter of convention whether we fix 7 to be a smooth measure 
on M, for example, or some power of a defining function times a smooth 
measure; any two such choices yield equivalent theories, but one does need 
to make an adjustment to the index set a below based on this choice. We 
shall follow the convention that 7 is instead a volume form for some fixed 
KQ metric, and hence is of the form x~ n ~ d times a smooth measure. This 
has the advantage that most of the operators in our later applications are 
self-adjoint. 

Basic facts about these spaces of pseudodifferential operators include the 
composition law 

(11) o C a" = (a w ,a' 01 ,0), 

which holds provided <7oi +cr' 10 > — 1 (this condition is needed to ensure that 
the integration defining the composition makes sense), and the existence of 
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a short exact symbol sequence 

(12) -► * ^q 1 '" tfgg S^( K ®T*M) -► 0. 

The other key facts we need concern the mapping properties of these oper- 
ators. To state these we first describe the appropriate function spaces. Fix 
any smooth K0 metric g on M. This determines the space L 2 (M; dV g ), as 
well as the basic Holder space A^q(M), which by definition is the closure of 
bounded 'W co functions with respect to the norm 

II II I / M \ U (P) — U (0)\ 

\\U \\o,a = SUp \u(jp) + SUp — • 

P eM dist 9 (»,g) a 

dist g (p,q)<l 

Next, for any s G N, set 

F| (M) = {u : Fi . . . V t u G L 2 (M; a%) : G r K0 V j < £, i < s}. 

Ajg(M) = {« : . . . V<« € A^e : V)- G r K0 V j < t, I < s}. 
and finally, for any defining function x for dM and 5 6 K, 

x S Hl e (M) = { u = x 5 v:vG H^ e (M)}, 

x»A^(M) = {u = x»v:veA^(M)}. 

There is a somewhat loose relationship between certain of these weighted 
Sobolev and Holder spaces. This is based on the fact that x a G x s H§ Le (at 
least locally near x = 0) if and only if a > 5 + (re + d — l)/2; similarly, 
x a G x^A^q near x = if and only if a > za Because of this we say that 
x^H^q and x^Aj^g are commensurable when u = 5 + (n + d — l)/2. 

Proposition 9. Fix fieN, o"io,ctoi,<5 G R smc/i i/iai aoi + <5 > 0, aio > (5. 
Lei A G *^'@(M). Taen i/ie maps 

A:x s H^(M) — x^ (M) 
A . s *+(n+d-i)/2 A fc+M.«( M ) _^ a; *+("+ d - 1 )/ 2 A^(Af) 

are bounded. 

We sketch a few points in the proof. First note that since these spaces are 
defined relative to K6 derivatives, and since ifae ° ^ko c ^Ke^ we ma y 
immediately reduce to the case [i = 0. Furthermore, conjugating by x~ s in 
the first case and x ~^ + ^ n+d ~ 1 ^ 2 ^ in the second, and observing that x/x lifts 
to a conormal function on Mj|q which is smooth (and nonvanishing) up to 
-Bn reduces us further to the unweighted case. Finally, the boundedness of 
^^q(M) on L 2 (M; dV g ) may be deduced via Hormander's method of using 
the symbol calculus to find B G satisfying A* A + B*B = C 2 Id + R 
for some C > and R G and then using Cauchy-Schwarz to prove L 2 

boundedness for operators of order — oo. The boundedness of on Aj^q 
is also deduced in two steps. Decompose A into a sum A' + A" where the 
Schwartz kernel of A' vanishes to infinite order at B\q and Bqi and that of 
A" is smooth across the diagonal. The boundedness of A' on Holder spaces is 
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equivalent to the standard local boundedness of pseudodifferential operators 
(of order 0) on a neighbourhood in M n , cf. [27]. This argument is discussed 
in detail in [19] and [17] for the case K = R. 

Proposition 10. Suppose that A £ ^^(M), where a = (oio, o"oi, cm) 
indicates conormal order of vanishing at each of the three boundary faces 
Bio, Bqi, B\\. If /j, < and an > 0, and if o§\ + 5 > 0, cio > 5, then A is 
compact on x s L 2 and on x <5+ ^ n+a!_1 - ) / 2 Aj^Q. 

This follows directly from the Arzela-Ascoli theorem. 

The parametrix construction for fully elliptic operators. We finally 
apply the theory of K0 operators outlined above to prove that under cer- 
tain hypotheses, the linearized gauged Einstein operator Laplacian L 9 is 
an isomorphism on certain weighted L 2 and Holder spaces. We state the 
result in slightly greater generality for an arbitrary generalized Laplacian 
P = V*X7 + R G Diff| (M; E), associated to a KQ metric g. The symmetry 
of this operator with respect to the volume form dV g simplifies some of the 
numerology below, but all of these results have direct analogues for more 
general fully elliptic IK0 operators. 

Definition 11. The generalized Laplacian P 9 is fully elliptic if its K0 sym- 
bol is invertible as a section of < Tf°°( KQ T*M; Hom(i?)), and if, in addition, 
its normal operator N(P 9 ), which is identified via Proposition^with the cor- 
responding operator P 90 on the model hyperbolic space ¥LH m , is invertible 
as an unbounded operator on L 2 (KH m ;Eo). 

Before stating the main theorem of this subsection, let us explore the 
relationship of this full ellipticity condition with the indicial root structure 
of P 9 . The indicial roots of P 9 are the roots of the indicial polynomial for 
P 9 , and give the rates of vanishing of formal solutions of this operator. If £ is 
an indicial root, then there exists some <}>{y) such that P(x^4>(y)) = &{x^ +1 ). 
The indicial roots of P 9 and of its normal operator N(P 9 ) are the same. The 
indicial roots are arranged symmetrically around (n + d — l)/2 in C. The 
complement in M of the set of real parts of all indicial roots of P is a union 
of open intervals and half-lines, again symmetric around (n + d — l)/2. The 
significance of these intervals is as follows. First, if 8 + (n + d — l)/2 is the 
the real part of some indicial root, hence at the boundary of two contiguous 
intervals, then neither of the mappings 

(13) P:x S Hl e {M,E;dV g ) — x 5 L 2 (M,E;dV g ) 

(14) P :x &+ ^ n+d - l ^ 2 A^{M,E) — > x^+^^A^M, E) 

have closed range. This is straightforward to check from basic definitions. 
Significantly deeper is the 

Theorem 12. Let P 9 G Diff^ (Af, E) be a fully elliptic generalized Lapla- 
cian. Define 5o by the condition that (n + d — l)/2 ± <5o are the real parts 
of the indicial roots of P closest to (n + d — l)/2. If 5q > and \5\ < Sq, 
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then (13\) and ( fLfl ) are both Fredholm, and are isomorphisms if and only if 
the nullspace of P 9 on L 2 (M;dV g ) is trivial. 



Proof. The fact that these mappings are Fredholm when \S\ < 5q will fol- 
low immediately if we can establish the existence of a parametrix G G 
y-^ a >°(M;E) for P with the property that PG - I = GP - I = Q G 
' a '°°(M ; E). Here for convenience we set r = (n + d — l)/2 + 5q and 
o~ = (fiOifoi) = ( T i r )- The final index (here or oo) corresponds to 

This parametrix is constructed in stages. We first choose an element Go in 
the small calculus, i.e. G G *^|(M; so that PG = I-Qo, Qo G • 
This uses only the symbol calculus and the symbol ellipticity of P, and 
proceeds exactly as in the usual (local) elliptic parametrix construction. 

For the second step we seek a correction term G\ G chosen so that 

the remainder term Q\ = I — P(Gq + Gi) lies in ^jjq ' a ' , and in particular 
is compact. For this we must solve the normal problem N(P)N(G\) = 
N(Qo) G ^(KH m ;E). By the second part of the full ellipticity hypothesis 
there is a unique solution to this equation in L 2 (WLH m , E) and (since the right 
hand side is ^q ) it is a simple matter to check that the solution is conormal 
at the boundaries of the quarter-sphere fibres of the front face. Indeed, using 
the analysis from |2|. we obtain that N(Gq) G -stf a XS++ 1 ), where a = (r, r) 
gives the orders of conormal vanishing at the two boundaries of the quarter- 
sphere fibres. 

Using |£| < Jo, we see from Proposition [9] that G\ and Q\ are bounded 
between these weighted spaces; from Proposition [10] we obtain also that Q\ 
is compact. This already shows that P 9 is Fredholm. However, it is useful 
to refine this parametrix further. 

Using the composition formula for K0 pseudodifferential operators, we 
see that the iterated compositions of this error term with itself vanish to 
increasingly high order at the front face, specifically Q\ G ^■^b ' 7 '" 1 • ^ e can 
therefore take an asymptotic sum of the series R ~ Yl'jLi Qi as an element 
of ' <T ' 1 - Now multiply P(G + Gi) = I- Qi on the right by I + R. We 
see that G = (G + G 1 )(I + R) satisfies PG' = I - Q' where Q' G ■ 

This error term lies in the very residual space of smoothing operators 
with Schwartz kernels which are conormal on M 2 . One consequence is that 
elements of the nullspace of L in either of these function spaces are conormal 
and vanish like x T . Furthermore, since these very residual operators form a 
semi- ideal (on one of these weighted I? spaces, say), a standard argument 
(cf. [19]) shows that the true generalized inverse G of any of the maps be- 
tween weighted spaces, which a priori is only defined as a bounded operator, 
is actually an element of , I / K Q Cr ' ; the error term PG — I = GP — I = Q is 
the projector onto the nullspace and is still very residual. If the nullspace is 
trivial, then Q = and hence P is invertible. □ 



22 



OLIVIER BIQUARD AND RAFE MAZZEO 



We remark, but do not prove, that if 5 lies in any of the other open 
intervals or half-lines described above, then these maps have closed range but 
are not Fredholm since either the kernel or cokernel is infinite dimensional. 

With not much more effort, we can prove that the Schwartz kernel of G 
has a poly homogeneous expansion at all boundary faces of Mj|q (as well as 
a polyhomogeneous expansion along the lifted diagonal of this space which 
is smoothly extendible across the front face). This implies that G maps 
polyhomogeneous sections to polyhomogeneous sections, and also shows that 
any section k which satisfies Lk = (even just in a neighbourhood of infinity) 
must have a complete polyhomogeneous expansion there. 

The resolvent family. The invertibility of the linearized gauged Einstein 
operator L 9 on weighted Holder spaces, which is a direct consequence of 
Theorem [12] and Proposition [9J is the key ingredient in the deformation 
theory of AIKH Einstein metrics. For the analogous result on products of 
AIKH spaces, we shall use a spectral synthesis formula for the inverse of 
this operator which expresses it in terms of the resolvent families of the 
corresponding operators on each factor. In preparation for this, we now 
recall the basic theory of these resolvent families in terms of the K0 calculus 
and prove some estimates on their Schwartz kernels which are uniform in 
the spectral parameter. 

First let us recall some results from [2] § 1.4]. Define Sq = 5^ for the 
operator L 9 as in the statement of Theorem \12\ then 

, . So =(n + d-l)/2 for IK = R or C, 

^ ) S > (n + d - l)/2 for IK = M or O. 

This means that the interval of weights 5 for which L 9 is Fredholm is exactly 
(0, n + d — 1) in the real and complex cases, and is larger in the quaternionic 
and octonionic cases. 

By definition, the resolvent of L 9 is the family of 1? bounded operators 
{L 9 — A) -1 , which exists precisely when A ^ spec (L 9 ). We wish to recognize 
these operators as elements of Vl^@ (M), depending holomorphically on A 
in an appropriate sense. This will follow from Theorem [T2l which in turn 
requires the 

Lemma 13. The operator L 9 — A is fully elliptic (as a IK0 operator) if and 
only if A ^ [<5q,oo). 

Proof. First note that 

Ke a 2 (L 9 - A) = K0 a 2 (V*V) = |C| 2 Id, 

which is obviously invertible. In addition, N(L 9 — A) = L 90 — A, so we 
conclude that L 9 — A is fully elliptic if and only if A ^ spec (L 90 ). 

The indicial operator of L 90 is a second order matrix-valued ordinary 
differential operator, and the indicial roots correspond to solutions of the 
form x^kq, where kq is a constant symmetric two-tensor. By reducing to 
the various irreducible components in Sym 2 , we obtain them as the roots 
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of a finite number of quadratic polynomials £ 2 — (n + d— 1)£ + a, where 
a is a constant depending on dimension and the irreducible component of 
the decomposition. The roots from any one of these polynomials are (n + 
d — l)/2 + + d — l) 2 — 4a. For some ao we obtain the roots with real 

part closest to {n + d — l)/2, that is (n + d — l)/2 ± <5o. So we see that 
each a < ao- Now, the indicial roots of L 90 — A are the roots of the various 
polynomials £ 2 — {n + d— 1)£ + a + A, hence are equal to 

{n + d- l)/2± \\/{n + d - l) 2 - 4a - 4A. 

Define 

<5 (A) = min \v^J{n + d - l) 2 - 4a - 4A = Re \[$k - A. 

a 2 V 

By Theorem 1121 -^ 9o — A is at least Fredholm on L 2 provided 

<$o(A) > 0, 

or equivalently, if A ^ [5 2 ,oo). 

So far we have proved that if A is outside this half-line, then L 90 — A is at 
least Fredholm. This shows that spec {L 90 ) is the union of [5 2 , oo) and finite 
point spectrum of multiplicity in the half-line (— oo, Sq). However, this point 
spectrum must be empty, since otherwise, if L 90 cj) = Xcj), for some A < <5 2 and 
with G L 2 , then the subspace spanned by all translates of <f> by isometries 
of ¥LH m would be infinite dimensional, contradicting the fact that L 90 — A 
is Fredholm. This finishes the proof. □ 

The same reasoning leads to the 

Theorem 14. Let L 9 be the linearized gauged Einstein operator on the 
manifold M with K@ metric g. Then L 9 — A is Fredholm if and only if 
A ^ [Sq, oo). More precisely, 

spec{L 9 ) = [5loo)U{X t }t 1 , 

where A.; lies in {— oo,<5 2 ) and is an 1? eigenvalue of finite multiplicity. 

An AKH Einstein space (M, g) is nondegenerate if and only if is not in 
this point spectrum. 

We shall need to know slightly more about the dependence of the inverse 
on A. 

Proposition 15. Fix e > and define Q e C C \ [5q, oo) to consist of the set 
of all A for which 5o{X) > e. Let t £ = (n + d — l)/2 + e and a £ = (r e ,r e ,0). 
Then for each e > 0, the resolvent family 

Q £ 3X^ R{X) = {L 9 - A)" 1 G ^^{M; Sym 2 ( K0 T*M)) 

is meromorphic in the sense that the Schwartz kernels of these operators, 
as elements of a fixed space of distributions, depends meromorphically on A. 
The poles occur only at each Aj and these are all simple; the residues are the 
finite rank orthogonal projections onto the corresponding eigenspaces. 
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The proof is based on the fact that the model resolvent (L 90 — A)" 1 is 
itself holomorphic, which can be checked by direct ODE analysis, and the 
analytic Fredholm theorem. This is the direct generalization of [25] and [9j, 
cf. also P2]. 

We conclude this section by proving uniform estimates for the off-diagonal 
Schwartz kernel of this resolvent when A = i\i lies on the imaginary axis. 

Proposition 16. Let (M,g) and L 9 be as above. When the indicial 

root of L 9 — in which has the smallest strictly positive real part is equal to 

CO*) = — - + y/#o-W- 

In particular 7(//) '■= Re C(aO > n + d — 1 with equality if and only if jj, = 0. 
Furthermore, 

C(m) ~ v\fA/% (i =t as — > ±oo. 

Lei K(z, z' ,ifi) be the Schwartz kernel of this operator. Then there exists 
fj,o > swc/t i/iai /or any c > and e G (0, 1), «/ z,z' G M satisfy d(z, z') > 
c > and > /iOj we /iat>e 

where the constant C is independent of fi. A similar estimate holds for all 
b-derivatives of K. 

Proof. The calculation of the indicial root and the statements about its 
asymptotics are straightforward, based on the remarks in the proof of Lemma 
[131 As for the main assertion, first suppose that / G ^^(M), |/| < 1, and 
define u = u M = R{ifi)f . We claim that for any e > 0, there exists a constant 
C £ > such that for > fiQ, 

\u\ < C/x-V 1 -^). 

In particular, C is independent of (and in addition, sup \ u\ depends linearly 
on sup |/|). 

To prove this, first recall that since R(ifj,) G ^jje 7 ^' 7 (M) , it is im- 
mediate that \u\ < AjafW; the issue is to prove the uniformity in fj, of 
the constant A. Suppose this fails, i.e. suppose there exists a sequence of /i 
tending to infinity so that 

sup /i|u(;z)|ar (1 - <07( ' i) = A„ -» oo. 
This supremum is attained at a point G M, so if we define 
w(z) = (^y-^M,) u(s), 

then 

(16) Kz)| < (x/x^))^"^, 
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with equality at z = and 

(17) {n~ l L - i)w(z) = (^)- (1 - £) ^VV) /(*)■ 

We shall consider various cases depending on whether or not q^ remains in 
a compact set of M . 

Suppose first that x(q^) > c > 0. Then the right hand side of (fTTj) tends 
to zero uniformly. Let B be a geodesic ball of radius 1 centered at fix a 
trivialization of the bundle over B and suppose that z are Riemann normal 
coordinates in this ball. Set £ = y/Jiz, so that £ lies in a ball of radius y/Jl 
in R n+1 . In terms of these coordinates, 

i=i ^ 

The remainder term is a second order operator with coefficients which con- 
verge to uniformly on compact sets in these expanding balls in R n+1 . 
Using standard local elliptic theory, we can take a limit in (I17p . and obtain 
a function defined on the entire Euclidean space, such that 

(A ? - ijWoo = 0, | Woo | < 1, I Woo (0)1 = 1. 

However, no such function exists. To see this, take Fourier transform (for w 
as an element of and use that the full symbol |£| 2 — i is nowhere 

vanishing. Hence this case cannot occur. 

Now suppose that x(q IM ) — ► 0. Pass to a subsequence so that — ► q £ 9M, 
and then apply a sequence of parabolic dilations based at appropriate 
points converging to q and with strength x(q , /J ))~ 1 so that D^q^) = (1,0) 
in a fixed coordinate system z. We now proceed much as before. Let B be 
a unit geodesic ball centered at (1,0), and define £ = y/Jiz. The sequence of 
operators ^~ l D* fl L converge to as before. The bound onw = D*w now 
takes the form 

\w\ < e ^ 1_e ) 7 W/v^ 

with equality at the origin; here t = \og{x/x{q tl )), and we can assume this 
is the first coordinate £i in the £ system. We again pass to a limit. The 
limiting function u>oo satisfies |u)oo(0)| = 1, 

(A e -i)w oo = 0, < eC^/A 

To analyze whether this is possible, note that this exponential bound on 
Woo implies that its Fourier transform is well defined as an element of 5?' on 
the subspace {£ 6 C n+1 : Im£i = (1 - e)/\/2,^ G R, j > 1}. The symbol 
£ • £ — i is again invertible here, which precludes the existence of this limit; 
hence this case is also impossible. This proves that the function u = R(i(j,)f 
satisfies the stated bound uniformly in \i. 

An essentially identical argument proves that a similar bound holds re- 
gardless of the location of the support of /. In other words, suppose that 
SU PP (/) C Bi(p fl ) and sup |/| < 1. Then for the L 2 solution to (L — ifj,)u = 
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/, we have \u(q)\ < C/i~ 1 exp(— (1 — e)d(q,p^)). The only modification 
needed is that if there a sequence for which the constant increases with- 
out bound, and if the center of the support of / M tends to infinity, then we 
parabolically rescale so as to obtain a sequence of problems (L^—i^u^ = / M , 
where the rescaled operators converge to the limiting model operator for 
the K0 structure. The validity of the bound in this case follows by what we 
have done above. 

We have now proved that K(z,z',i{i) decays like e -( 1 ~ £ ) d ( z < z for 
d(z, z') > c > in a weak sense. More precisely, let °k be any neigh- 
bourhood in M x M with compact closure which does not intersect the 
diagonal; then if a > 2n + 2, the H~ a norm of the restriction of Schwartz 
kernel to satisfies this bound. Using that (L z — i\x)K = [L z i — ifj,)K = 
away from the diagonal, we can estimate any norm of K in ^ at the 
cost of introducing an extra factor /i fc+cr into the estimate. This in turn may 
be absorbed into the exponential by decreasing the factor e slightly. This 
completes the proof of the bound, and indeed also of bounds with respect 
to any c € k norm in the interior. In fact, it gives slightly more, namely that 
this bound holds even after applying any sequence of K0 vector fields to K 
on the left and right; this is because ICO derivatives are controlled by powers 
of L, which as above are equivalent to powers of \i. 

To finish, we also need to check the conormal bounds, i.e. that the same 
estimates remain true if we apply any sequence of b vector fields to K on 
the left and right. For this we point out the following facts: first, since 
K(z' ,z,ifi) = K(z,z', —ifj,)*, we need only the case where all b derivatives 
are applied to the left - z - factor; next, K itself is conormal, so these 
6-derivatives behave well locally uniformly in /x, i.e. it is only the large /x 
behaviour that might be problematic; finally, we can repeat the same proof 
as for K itself, using at the final step to convert the weak bounds to strong 
ones that if V is any 6-vector field, then [L, V] is a ICO-operator of order 2, 
hence is bounded by multiplication by /x □ 

Remark 17. We have stated the results on the resolvent family for the 
linearized Einstein operator, but the results remain true for any geometric 
Laplacian, provided we choose <5o as in Theorem [T2j For example, on an 
asymptotically quaternionic hyperbolic space, Theorem Q3] gives the spec- 
trum of the Hodge Laplacian acting on differential forms (except when the 
degree equals half the dimension, then there is a zero eigenvalue of infinite 
multiplicity) . 

4. Einstein deformation theory 

We now present some basic facts about the (Bianchi gauged) Einstein 
operator and its linearization. Using results from the last section, we review 
how this yields the deformation theory for the rank one hyperbolic spaces 
in the class of AICH Einstein spaces. This is contained in [2] for all K, see 
also [H] and [15] for the result when K = M, so the only novelty here is 
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showing how this follows immediately through the use of the K0 calculus. 
These same arguments are used again in § [6] for a coupled generalization of 
these same equations, and in the product case in § 

4.1. The Einstein equation and the Bianchi gauge. The Einstein 
equation Ric 9 + Xg = is not elliptic because of its diffeomorphism invari- 
ance. Amongst many viable gauge choices, the Bianchi gauge introduced in 
[2] is particularly convenient. Define the map from symmetric 2-tensors to 
1-forms, relative to the fixed background metric g, 

(18) k i — > B 9 (k) = 5 9 k + ^dtr 9 k. 

Note that B 9 (g) = 0, so the subspace of metrics g near to g which are in 
Bianchi gauge (i.e. so that B 9 (g) = 0), is identified with the set of tensors k 
near such that B 9 {k) = 0. The system 

(19) Ric(5) + A5 = 0, B 9 (g)=0, 

which is elliptic in the sense of Agmon-Douglis-Nirenberg, can be rolled up 
into the single elliptic equation 

(20) N 9 (k) := Ric(<7 + k) + X(g + k) + (5 9+k )*B 9 {k) = 0. 
As proved in [21 chapter 1], 

Proposition 18. Suppose that N 9 (k) = 0, and in addition that \B 9 (k)\ — > 
and the Ricci curvature of g + k is nonnegative and strictly negative some- 
where. Then g + k satisfies $19]) . i.e. is Einstein and in Bianchi gauge. 

The proof follows from the Weitzenbock formula 

B 9+k N 9 {k) = 5 9+k {5 9+k )*B 9 {k) = ((V 9+k )*V 9+k -Ric 9+k )B 9 (k) 

and the Bochner technique. 

From the same Weitzenbock formula, the converse follows almost imme- 
diately, that is, any Einstein metric g close to g can be put in the Bianchi 
gauge to satisfy the system (fT9|) . More precisely, let Diff denote the set of 
all diffeomorphisms on the AKH space M which are close to the identity 
and exponentials of vector fields X G x^AJ^q(M), ^ the set of all metrics 
g = g + k with k £ x^A^q(M), and 5? C ^ the set of metrics g which 
satisfy B 9 (g) = 0. Then one has the following slice statement [21 chapter 1]: 

Proposition 19. //Ric 9 < 0, then the natural map 

Diff x y — ► J(, (g, $) i — ► (g) 
is a local homeomorphism. 

An advantage of this gauge is that the linearization takes the simple form 

(21) L 9 k := 2DN 9 \ (k) = V*V k - 2Rk + Ric o k + k o Ric + 2 A k; 
here 

o 

(yRhx^jj^j — Ripjg , (R,1C O K/^jij — R,1C^ Kjpj j (a£ O R^C^)^ — R-lCpj, 



28 



OLIVIER BIQUARD AND RAFE MAZZEO 



and all curvatures and covariant derivatives are computed relative to g. Note 
in particular that if Ric 5 = —Xg, then 

(22) L 9 = V*V - 2R. 

4.2. Deformation theory for AKH Einstein spaces. We now review 
the basic deformation theory for AKH Einstein spaces, proved originally in 
[llj when K = R, and in [2] in the other two cases. We do not discuss the 
more subtle aspects of this deformation theory, but restrict attention to the 
simpler case of perturbations of nondegenerate AKH Einstein metrics. 

Definition 20. An A¥LH Einstein metric g is said to be nondegenerate if 
the L 2 nullspace of the linearized Bianchi- gauged Einstein operator L 9 is 
trivial. 

Recall that g is nondegenerate if for any v > —5^, the nullspace of L 9 on 
x u+ ^ n+d ~ 1 ^ 2 A^Q(M,S 2 (T*M)) is trivial. This follows from the regularity 
theorem stating that if L 9 k = and |k| < Cx v for some v' > (n + d — 1)/2 — 
<5o, then k 6 ^^^ d 1 ^ 2+<5 ° ) and in particular k £ L 2 . 

The significance of this nondegeneracy condition is contained in the 

Proposition 21 ([TT], [2], [15]). Let g be a nondegenerate AM>H Einstein 
metric with "r^ 00 conformal infinity c. Then every "if 00 conformal infinity 
datum c' sufficiently close to c in the c € 2,a topology is the conformal infinity 
of an AWLH Einstein metric g' such that g' — g c i G x^A 2 ^ for some v > 0. 
(Here g c i is the A¥LH metric with conformal infinity data from Definition 
2 in §2.) This metric g' is unique amongst AWLH metrics with the specified 
conformal infinity and such that \ \g' — 5 c '||2,a,z/ is small. 

Proof. First define an extension operator which associates to the conformal 
infinity c' an AKH metric g c i . If c = ( [7] , rj) , then we choose a product 
decomposition (0, e) x X of a collar neighbourhood of dM and a radial 

coordinate x so that g has the form dx x t^ + ^¥ + & 5 with k 6 xA^q. (We can 
take the weight v to equal 1 when c is smooth, but could also use any smaller 
positive value.) Fixing a cutoff function x( x ) which equals 1 for x < e/3 
and vanishes for x > 2e/3, then for any conformal infinity c' = ([7'],??')) se t 

, \ [dx 2 + i/A (r/) 2 \ 

ge = (l - x(x))g + x(x) + Tk + k J- 

(In the quaternionic case, the function spaces aAj^g vary with c' so one has 
also to choose k varying continuously with c'). It is straightforward that 

\\N 9 (g,)\\o, a ,u < C (|| 7 ' - 7 || 2)Q + \\r/ - r,\\ 2ia ) . 
Writing N 9 (g c i) = f c i, then Taylor expansion gives 
N 9 (g c ,+k) = f c , +L c 'k + Q(c',k), 

where the second term on the right is the linearized Bianchi-gauged Einstein 
operator at g c /. The nondegeneracy of L 9 implies that L c is also invertible 



A NONLINEAR POISSON TRANSFORM 



29 



if c' is close enough to c; we denote its inverse by G c . When c' is this 
operator is a IKO pseudodifferential operator of order —2, and for < v < 1, 

G c ' : x v A^q(M,S 2 (T*M)) — ► ^A^(M,5 2 (T*M)) 

is bounded. Furthermore, the norm of this operator is bounded indepen- 
dently of c' in a neighbourhood of c. When K = H, the function spaces vary 
with 7/ (i.e. with the distribution $>'). 
Now write the equation to be solved as 

k = -G e [ft + Q(c',k)). 

The right side defines a contraction mapping when c' is sufficiently close to 
c, and from this we immediately obtain a unique solution k. The metric 
g c / + k is an AKH metric which solves the gauged Einstein equation. By 
Proposition 1181 it is in fact an AKH Einstein metric in Bianchi gauge. □ 

Remark 22. A precise statement about the regularity of g and g' near dX 
has been omitted, and indeed this is a subtle issue. There is a substantial 
difference between understanding the dependence of the asymptotic regu- 
larity for an arbitrary AK.H Einstein metric on that of its conformal infinity 
data, and the same question for such a metric obtained by perturbation from 
one which is a priori known to be polyhomogeneous. The reason is that in 
the former case one needs to deal explicitly with the gauge conditions, while 
in the latter, the gauge choice is part of the setup and the perturbation term 
automatically satisfies a IK© elliptic equation. Here is a summary of what 
is known. 

For the nonperturbative case, when K = M and the conformal infinity 
data c of the ARH Einstein metric g is < ^°°, then g is polyhomogeneous, 
and in fact, in even dimensions has a smooth conformal compactification, 
[5], [13] . The corresponding result has not been proved in the other cases 
(except when IfC = C and g is Kahler-Einstein, [E]), but is surely true by 
essentially the same method as in [5]. 

As for the perturbative case, a simple adaptation of the argument in 
|20j . which depends only on the commutation properties of with b 

vector fields on X, proves that if the perturbed conformal infinity data c' 
is < ^°°, then the solution k, and hence the metric g' , constructed above is 
polyhomogeneous. If K = R and c' is only C 2,a , then [15] proves that if 
g is sufficiently regular, then g' has a c ia 2,a conformal compactification, but 
this is not a sharp statement since the function space ' a (dX) is not well 
adapted to these types of degenerate problems. 

These issues will not be emphasized here, and we shall tacitly assume the 
polyhomogeneous regularity of AKH metrics with < &°° conformal infinities. 
To ease the reader's conscience, however, since we will not supply the full 
proof of that fact, note that this issue only arises in § d and one can easily 
adapt the arguments to accommodate metrics with lower regularity, as we 
discuss briefly there. 
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5. Asymptotically product hyperbolic metrics and their 

conformal infinities 

Let M" l+1 = Gi/Ki be an AIQH space, i = 1, 2. The boundary at infinity, 
S ni = Ki/Hi, is equipped with the standard IEQ-contact distribution 
which has a conformal Hi structure inducing a compatible conformal class 
[ji]. The product hyperbolic space M = M\ x M 2 is a (reducible) rank two 
symmetric space with Furstenberg boundary X = S ni x S U2 . Let Qi be the 
standard metric on each factor, so that Ric 9i + \gi = 0. Then 

g = A151 + X 2 g 2 

is Einstein with 

Ri c f + g = 0. 

We begin this section by describing a class of boundary structures on 
X, called (G\ x G^-conformal structures, which constitute the conformal 
infinity data for the class of Einstein metrics we eventually construct. The 
problem of extending one of these boundary structures to a metric on M 
which is asymptotically Einstein in an appropriately strong sense is far from 
immediate. The 'obvious' extension has Einstein tensor vanishing in some 
sector near X, but not uniformly near infinity. The main goal of this section 
is to explore the geometry of asymptotically product hyperbolic metrics 
in order to find the correct compatibility conditions for metrics which are 
asymptotically Einstein in this stronger sense. Their construction is carried 
out in the next section. 

5.1. (Gi x G2)-conformal structures. Using the notation above, we make 
the 

Definition 23. A (G\ x G^-conformal structure on X consists of a pair 
of distributions, each equipped with a conformal class of metrics, (^, [7^]), 
i = 1,2, such that 

(1) the distributions &i = @i + are integrable; 

(2) ^9^2 = TX; 

(3) the pair (S^i,^/ 'S>i) with induced bracket 

[ , ] : % x % ^ 

is isomorphic to the graded Ki-Heisenberg algebra; 

(4) Q>i is equipped with a conformal Hi structure, compatible with the 
bracket, and inducing the conformal metric [7*]; 

(5) [^1,^2] C ^i + ^2- 

The basic example, of course, is a product structure on X = X\ x X2: 
here fiC#i = TX 1 © {0} and Q> 2 C ^2 = {0} © TX 2 , and each {%, [ 7i ]) 
is the pullback of a Gj-conformal structure from X^. Our main focus in this 
paper is with perturbations of these product structures. As we now indicate, 
there is substantial rigidity in the deformation theory, and nearby structures 
retain many vestiges of the product case. 



A NONLINEAR POISSON TRANSFORM 



31 



Lemma 24. For i = 1,2, let Xi be a compact simply-connected manifold 
with Gi-conformal structure (f^,[7i]). Then any small deformation of the 
product (G\ x G2) -conformal structure on X = X\ x X2 has the following 
properties: 

(1) the pairs of distributions (3>i,3>2) and (#1,^2) remain of product 
type; 

(2) if Kj = R or C, i/ien /or that i, Q>i remains fixed (up to a global 
diffeomorphism), but the deformation of [fi] may depend on both 
factors in X\ x X^; 

(3) if Kj = H, i/ien the distribution S>i varies amongst quaternionic 
contact structures, but since the conformal class [73] is determined 
by Q>i , any deformation of [7^ depends only on the factor Xi ; 

(4) if Kj = O, then both 3>i and [7^] remain fixed in the deformation 
(and in fact Xi = S ni and the structure is standard). 

In particular, the distribution ^ can change (modulo diffeomorphisms) only 
in the quaternionic case, and the conformal metric [ji] may depend on both 
factors X\ and X2 only if Xi is real or complex. 

Proof. At the initial product structure, the leaves of the foliation corre- 
sponding to ,j^i are just the parallel copies of Xi. After a small deformation, 
the leaves are covering spaces for Xi, and since both X\ and X2 are simply- 
connected, these leaves must remain diffeomorphic to Xi. In particular, 
the perturbed distributions ^ are still equal to the tangent bundles of the 
respective factors. 

Now observe that as an immediate consequence of conditions (3) and (5), 
[&i,@2] = {%. + [^1,^1], ^2] C + ^ 2 

(by the Jacobi identity), and hence S>2 is invariant along the leaves of the 
foliation corresponding to #1; similarly S>\ is invariant along the leaves of 
the foliation for #2- 

When Kj = R, ^ remains equal to the tangent bundle TXi, while if 
HQ = C, then by Darboux's lemma, we may still assume that remains 
fixed in the deformation. 

In other words, in this deformation theory, we may as well assume that the 
distributions % and ^ remain of product type. The remaining assertions 
follow directly from this. □ 

When X\ or X2 are not simply connected, we shall impose this product 
structure as a separate hypothesis: 

Definition 25. A deformation of a product (G\ x G2) -conformal structure 
on X = X\ x X2 is called globally integrable if (modulo diffeomorphism) 
the foliations &\ and ^2 remain the tangent spaces of the two factors of 
X = X 1 x X 2 . 

It is possible to define (G\ x G^-conformal structures on any closed man- 
ifold X of the appropriate dimension. Looking ahead to the main goals of 
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this paper, one could then try to extend this to an asymptotically Einstein 
metric on some manifold M with two boundary hypersurfaces F\ and F 2 
and X as its corner of codimension 2. However, we have already noted that 
this extension problem is not at all easy; in fact, the main difficulty seems 
to be the extension from X to the boundary faces F{. For this, it appears to 
be almost necessary that X and M be products, and that the metrics and 
boundary structures are globally rather similar to the ones considered here. 
Thus, in all that follows, we shall assume for simplicity that X, M and the 
distributions Q>i and ^ are products. 

5.2. Asymptotically product hyperbolic metrics. It is always possible 
to construct a complete metric on the interior of M which is 'weakly' product 
hyperbolic and with any given (Gi x G2)-conformal structure on X as its 
prescribed conformal infinity. In fact, we can just write down a formula 
which directly generalizes (JH): let x, be defining functions for the boundary 
hypersurfaces X{ C Mj, and choose ImKj-valued 1-forms rji defining % 
and compatible metrics 7, representing the given conformal classes. In a 
neighbourhood of X of the form (0, e) Xl x (0, e) X2 x X, set 

dx^ r *f\ dx^ ^2 ^2 

(23) 9jim>l2 , m = Xi{^ + ^ + ^4) + + 4^ + 4^)- 

Slightly more generally, consider 

(24) g = 571,^1,72^2 + k, ke {xix 2 ) u C 2 ' a 

for some v > 0. (The norms and covariant derivatives are with respect to the 
metric (I23p .) We call any such g a weakly asymptotically product hyperbolic 
metric, and say that {S>\, S>2-, [71], [72]) is its conformal infinity. Just as in the 
rank one setting, g determines this conformal infinity. Conversely, replacing 
7i and rji in (|24|) by any other conformal representatives 

7i = fili, m = fitji, 
where the f\ are strictly positive smooth functions on X, yields a new metric 

~ _ x ( dx l , hn , fhl\ , x ( dx l , hi2 fhl\ 

We claim that up to a diffeomorphism g is asymptotically equivalent to 
g. Indeed, if &*Xi = and $\ x = id, then 

dx i dx <i 



dxi dx,i dfi 
— = wr, and 

X% Xi ZJi 



Xi X{ 



+ X 2 



and hence 



\$*g-g\ g = &(xi + x 2 ). 

A metric g as in (|23p is sometimes also called weakly asymptotically Ein- 
stein, because of the 
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Lemma 26. The metric (23\) satisfies the estimate 

\Ric 9 + g\ g = @{x x +x 2 ). 



The proof is deferred to the next subsection, where the formalism for the 
necessary calculations is developed. 

As indicated earlier, we also define a narrower class of strongly asymp- 
totically product hyperbolic (or strongly asymptotically Einstein) metrics, 
for which the Einstein tensor decays uniformly near the entire boundary of 
M. The goal in the next few subsections is to find the equations which the 
limiting values of g\ and 72, 772 must satisfy at x 2 = and the limiting values 
of 52 and 71,771 must satisfy at x\ = 0, in order that g lie in this smaller 
class. 

Later in the paper, in § 17.21 we shall also define a class of 'near product 
hyperbolic' metrics; these will be defined by slightly different conditions, but 
we show there that any strongly asymptotically product hyperbolic metric 
is of near product type. 

This profusion of similar names is indicative of the fact that for metrics 
which are modelled by symmetric spaces of rank greater than one, it is 
by no means clear what the precise conditions are under which a metric 
should really be considered 'asymptotically symmetric'; each of the classes 
of metrics above has some claim to this moniker in the product hyperbolic 
setting. 

5.3. Asymptotic curvature calculations. We now calculate the asymp- 
totics of the Ricci curvature for a weakly asymptotically product hyperbolic 
metric (|24|) on M\ x M 2 . In the course of this the proof of Lemma [26] 
will emerge, as well as motivation for the extra conditions imposed on g to 
warrant the name strongly asymptotically Einstein. 

As before, assume that M, X and the distributions 9>i C are all of 
product type. The main calculations are local near the boundary faces; to 
be definite we work in the region where x 2 — > 0, and write g in the form 



where the two terms are metrics along horizontal and vertical slices, M\ x 
{p 2 } and {p\} x M 2 , respectively. We assume that 

• the x 2 dependence is only what is written explicitly; in other words, 
gi , rj 2 and 72 are defined and smooth on Mi x X 2 and are independent 
of x 2 ; 

• rj2 is the pullback of a contact form from X 2 , hence is independent 



• 72 is a family of metrics on *2s 2 compatible with 772 (and hence gives 
a GVconformal structure on each slice {pi} x X 2 ). 



(25) 




of Mi; 
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The precise form of g\ is not so important for the moment, but in order to 
maintain consistency with (|24ft , we also impose that in analogous coordinates 
near the boundary of Mi, g\ ~ \\{dx\jx\ + 7i/4xf + "qf/Axf) as x\ —>■ 0, 
and that 72 converges to a representative of the specified conformal class [72] 
as x\ — ► 0. However, these last conditions do not enter into the immediate 
considerations. 

Well-known formulae due to O'Neill, cf. Proposition 9.36 in [T), express the 
Ricci curvature of a Riemannian submersion in terms of the Ricci curvatures 
of the base and fibres and two additional tensors: the second fundamental 
form T of the fibres and another tensor which measures the deviation of the 
horizontal subspaces from being integrable. To adapt this to our setting, 
we regard M as a fibration Mi x M2 — * M\. The two factors are orthogo- 
nal, and the horizontal subspaces are integrable (with leaves the Mi slices, 
i.e. the submanifolds Mi x {^2})) but this is still not quite a Riemannian 
submersion because g\ depends also on M2. In the curvature computations 
below, however, it behaves asymptotically as xi — > like a Riemannian sub- 
mersion: the negative powers of X2 in all terms in g<i add an extra X2 factor 
to all derivatives in the M2 directions. 

We continue by defining the various quantities which appear in the O'Neill 
formula?, and developing some of their properties. The first is the second 
fundamental form for 72. This is the section T of T*M\ <S> Sym 2 (^|) defined 
by 

(26) <T(&,C2Ui> = -^G%72)(6,C 2 ) **T = -\d Mx l2- 

(Here and later, a subscript 1 or 2 of a vector indicates the factor to which 
it is tangent.) For either of these expressions we regard 72 as a section of 
Sym 2 (f2?2*), which in turn is a trivial bundle over each Mi slice. For each 
(2 £ f 2 we also set 

(27) T 6 G End(^ 2 , TM\ ) , T| 2 G End(TMi ,@ 2 ), 

where the metrics g\ and 72 are used to dualize. 

There is still a freedom in the choice of the representative 72 G [72], but 
we now fix the normalization that the volume form dV^ 2 is constant in the 
Mi directions. Consequently, the mean curvature vector vanishes: 

Xr 72 T _ Q 

The trivial connection d 1 on &2 is not compatible with the metric, but 
to find one which is it suffices to add the second map in (|27j) : thus 

(28) V = d Ah + T, (i.e. V 6 6 = £^6 + T£&) 

defines a unitary connection on S>2 over Ml. The divergence of T with 
respect to this connection is the bilinear form on 3>2, 



(5 Ml T)(6,C 2 ) = -^((V eQ T)(6,C 2 ),e Q ), 
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where {e Q } is an orthonormal frame for TM\. The final ingredient we need 
is the bilinear form Q on TM\ defined by contracting the product of T with 
itself in the Sym 2 ^!) component with respect to 72, 

(29) Q(ti,Vi) ■= (T^,T m } 12 . 

Before proceeding, we derive the crucial first-order properties of T. 

Lemma 27. Let dy 1 denote the exterior derivative on M\ coupled to the 
connection V on S>2- Then 



d^ h T = 0. 



Proof. By (|2SJ), d^T = d Ml T + T A T, where in the last term we regard 
the two factors as elements of T*M 1 <8> End(^ 2 ) and T*M X <g> Sym 2 (^ 2 ), 
respectively. The second expression for T in (|26D gives d Ml T = 0, hence it 
suffices to prove 

T AT = 0. 

This identity in turn is a direct consequence of the symmetry of the action 

(if • q)(x, y) = q(ux, y) + q(x, uy) = 72 ((uu + vu)x, y), 

of a symmetric endomorphism u on a quadratic form q(x,y) = 72 (ux, y) in 
the pair (u, v). □ 



Lemma 28. On each slice Mi x {^2} there is a Bianchi identity of the form 

1 
2 



(S Ml Q + -d Tr Q) = (5 Ml T,T*0, £ £ TMi . 



Proof. Choose an orthonormal frame {e a } for TMj and extend £ to a vector 
field on Mi which is parallel with respect to V at some point z\. Then, 
calculating at Z\, 



(5 Ml Q)(0 = -^V ea Q(e a ,0 

a 

= -^{(y ea T*)e a ,T*0 + (T*e a , (V ea T*)£) 

a 

= {5 M 1TiT *Q_J2(T*e a ,(ytT*)(e a )) 

a 

= (5 Ah T,T*0-\dTrQ(Z). 

The second equality uses that dy h T = 0. □ 

There is also a second fundamental form Ij and corresponding mean cur- 
vature vector iVj = Tr 9i Ij for each Mi slice, i = 1, 2. Note that ./V2 is different 
from Tr 72 T (which we are assuming is equal to 0), since in the latter one 
only takes the trace in the £^2 directions. 

We can now state an exact formula for the Ricci curvature. 
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Lemma 29. Let (M = M\ x M%, g = <7i + 92) a metric on M keeping the 
factors Mi and M2 orthogonal. Let Ij be the second fundamental form of Mi 
and Ni the mean curvature vectors. Then the Ricci tensor of g is given by: 
• 

Ric 9 (6,Ci) =Ric 9l (6,Ci) - (<5 M2 Ii)(£i,Ci) - <Ii(£i,Cl),iVl) 

+ (5 Mi )*iv 2 (ei,Ci)-(i^i,iKi) 5 

wii/i an analogous expression for the restriction of Ric 9 to TM2. 
//ere 5* is the symmetrization of the covariant derivative. 

+(5 M2 i 2 (e 2 ),6) + (v 6 iv 2 ,ei). 

/fere <5 Mi Ij is i/ie divergence o/Ij regarded as a symmetric 2-tensor 
along Mj. 

The derivations of these two formulae are left to the reader. 
As a first application, we have the 

Proof of Lemma\M We will be applying ([3D]) and ((3TJ) with M { = (0, 1)^ x 
Xi. The second fundamental form of the slices M\ x {p 2 } is 

<Ii.6> = "^651 = -y (^2^671 + ^6^)- 
If £2 is a un h vector in TX 2 , then 

^2^671 + ^(.A = ^2), 

since ^1 does not depend on M2. Hence on each slice M\ x {£2}, Ii = &{x 2 ), 
and the same is true for all its derivatives. On these same slices one also has 

Ric 91 = -Ai<7! + G(x x ). 

Analogously, on the slices {pi} x M2, we have 

I 2 = 0(x x ), and Ric 92 = -A232 + 0{x 2 ). 

Inserting these in ([50]) and (f3Tj) gives Ric 9 = — g + G[x\ + x 2 ), as desired. □ 

The main result of this subsection is the 
Lemma 30. Let g be defined by $25\) . and suppose that 
( S Ml T = 

{ Ric 9l (6,Ci) + Ai(£i,Ci} 9l = <T*£i,T*Ci> 72 

for all vectors £1,(1 S TM%. Then 

Ric 9 = -g+ ff(x 2 ). 

In particular, if g is as in (2~4\ ) and the equations [32]) are satisfied at both 
boundary hyper surf aces, then g is strongly asymptotically Einstein. 



(32) 
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Remark 31. Note that the right hand side of the Bianchi identity in Lemma 
[28] vanishes when the first equation in (|32p is satisfied. 

Proof. We apply (p0|) and ([31]) as follows. First, just as before, Ii = G{x 2 \ 
and hence N\ = 0(x 2 ) too. On the other hand, the normalization on the 
volume form implies that A2 = 0, thus 

Ric 9 (Ci,Ci) = Ric 9l (6,Ci) - O^i^Ci) + 0(x 2 ), 
Ric 9 (6,6) = ^(x 2 ), 

Ric 9 (6,C2) = Ric 92 fe,C 2 ) - (5 A/l I 2 )fe,C 2 ) + G(x 2 ). 

From the formula (|25p and the fact that r\ 2 is constant along M\ slices, we 
get 



and since the connection (|28p on $> 2 along M\ is exactly the one induced by 
the Levi-Civita connection of g, the result follows. □ 

To conclude the section, observe that the formulae in the lemma corre- 
spond exactly to the standard formulas obtained for a Riemannian submer- 
sion with integrable horizontal distribution [H proposition 9.36], as expected 
from our claim that the asymptotic behaviour when x 2 — > is that of a Rie- 
mannian submersion. 



6. Extending the approximate solution to the codimension one 

boundary faces 

Let (M = Mi x M 2 ,g = + g 2 ) be a product of AKH Einstein metrics, 
with conformal infinity c° = (3^,0$, [7J], [7^]) on X = X\ x X 2 . As in 
the last section, we consider deformations c of c° (assumed to be globally 
integrable in case either X\ or X 2 is not simply connected). According to 
Lemma [24[ the pair of distributions ^1,^2 remains of product type on X, 
and we then extend these by pullback to a pair of transverse distributions 
of product type on all of M. Choose metrics 7$ representing each of the 
conformal classes [7$]. In the real or complex case, these may depend on 
both factors of X, but we maintain the normalization so that, still just over 
X, dV 71 is independent of X 2 , and similarly dV 12 is independent of X±. 
Based on the calculations of § 15.31 we now address the problem of how to 
extend (71, 72) over the faces Mi x X 2 and X\ x M 2 to obtain a metric which 
is strongly asymptotically Einstein. 

6.1. Extension along boundary faces. We focus on the extension of 71 
over Mi x X 2 , since the other case is treated exactly the same. 

Let us restate the problem more carefully. On the face Mi x X 2 , we seek 
metrics g\ on TM\ and 72 on 2> 2 which solve the system (|32|) . The solutions 
are constrained by the requirements that dV 72 is independent of Mi, and 
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that (gi ,72) is asymptotic at X = dM\ x X2 to the given (G\ x G2) -conformal 
structure in the sense that 



as x\ — > (with an error term 0{x\) for some ^ > 0). 

Lemma 32. Suppose that K2 = HI or O. T/ien (CHjj reduces to a single 
uncoupled equation on M\ which is simply the usual Einstein equation. 

Proof. Under this hypothesis, 72 is independent of the X\ factor. Hence 
T = and the first equation in (|32|) is satisfied. The second equation reduces 
to the uncoupled Einstein equation on M\. By Proposition [21] we can extend 
the conformal class [71] on X\ to an AK4H Einstein metric g\ on Mi; note 
that this is actually done parametrically, depending on q<i £ X2. □ 

When M2 is real or complex, (j32j) cannot be reduced in this way, but 
fortunately, solutions can still be obtained near to the standard one by per- 
turbation methods. 

We can now state and prove the main result of this section. 

Theorem 33. Suppose M = M\ x M2 is a product of AKH Einstein spaces 
such that the L 2 nullspace for the linearized gauged Einstein operator on Mi 
vanishes. Then, for any small globally integrable perturbation of the product 
(G± x G2)-conformal structure on X, the system has a global solution 
(91 5 72) on the face Mi x X2 with the prescribed asymptotic behaviour at X, 
more precisely on each slice Mi x {q2}, 



with smooth dependence with respect to 02- 

Proof. This proof is similar to that for Proposition [2T1 Consider the slice 
Mi x {q2}, and begin with the conformal structures 71, 72 on Xi x X2, with 
dV~* 2 independent of X\. Fix a smooth extension map assigning to 71 a 
metric gi on Mi with 



as xi — > 0. Here A 2 '" is the geometric Holder space on Mi x {(72}. The 
weight v is positive; we can fix v = 1, but any smaller value is possible. 
Recall that we already have 



Ric^i + gi = 0{xi). 
As before, extend 72 by pullback on Mi, so that the corresponding second 




72 



72, 




72 - 72 G < ^ 2 ' 1 



5i = Ai(^ + -^ + ^)+A;, k€(xix 2 rX 
x\ 4x| 4x| 



fundamental form T 
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When K = R, we consider perturbations (ft of 72 which fix dV 12 , so the 
tangent space consists of trace-free symmetric two-tensors, i.e. sections of 

y = Sjml(TX 2 ). 

When K = C, 72(v) = drj 2 (-,J-), where J is an almost complex struc- 
ture; the perturbation (ft must also remain compatible with dr\ 2 on 3t 2 i or 
equivalently is a deformation of J, so the tangent space consists of trace-free 
J-skew-Hermitian symmetric two-tensors, i.e. sections of 

^ = Sym (^ 2 ). 

Consider elements 

h e A^' a (Mi x {g 2 },Sym 2 (rMi)), (ft G A 2 L ,' a (M 1 x {q 2 },f), 

and assume that both h and (ft have sufficiently small norm. Suppose that 
(gx + h, 72 + </>) is a solution of the system ([32]) along Mi. 

Denote by and the second fundamental form and corresponding 
quadratic form defined by 72 + (ft. To break the diffeomorphism invariance 
of the equation, we add the Bianchi gauge condition 

B 9l (h) = (5 91 + ^dTi 9l )h = 0. 
Thus we consider the system 

&*™(h,<t>) = (Ric 9l+h - X( gi + h)-Q <t> + {5 9l+h )*B 9l {h),5 9l+h ^). 
By Lemma [28l any solution of this equation must also satisfy 

B 9 ^ +h (5 9l+h )*B 9l (h) = 0. 

By the same argument as in the uncoupled case, cf. Proposition [TBI we 
conclude that B 9l h = 0. Hence a solution of <I> 9l ' 72 (/i, (ft) = is also a 
solution of the original system f|32|) . 

As in the proof of proposition [2T1 it suffices to check that 

$91,72 . A£ a (Mi x {q 2 }, Sym 2 (TMi) © S>) 

— A°' a (M! x fe}, Sym 2 (TMi) © ^) 

is a "r^ 1 mapping of Banach spaces, for (h, (ft) of sufficiently small norm, and 
furthermore, that its linearization 

D^\ m : A^Mi x {q 2 }, Sym 2 (TMi) © y) 

— A°' Q (Mi x to}, Sym 2 (TMi) © J?") 

is an isomorphism at the product metric. Since the linearization of the 

o 

Bianchi-gauged Einstein equation is V*V — 2R, and T = at the product 
metric, this linearization decouples as 

D$| (00) (h,<ft) = ((V*V - 2R)h, V*V<ft) 
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for the metric g®. By the hypothesis on the vanishing of the L 2 nullspace 
for (Mi , g\ ) , the first component is an isomorphism. The second component 
V* V is the rough Laplacian, and this is an isomorphism for weights v G 
(0,n + d-l). 

The last statement comes from the smooth dependence of the solution 
constructed by the inverse function theorem with respect to the parameter 
1i- □ 

Corollary 34. With the same hypotheses as in Theorem \33\ let (91,72) be 
the solution of the system liSty) on the face M\ x X<i- Then the metric 

h= +A ( ^ I ^ 2 I V<2 
\ x\ Ax?, 4x2, 

defined in some neighbourhood of the face M\ x X2 where X2 <C 1 on M\ x M2 
satisfies: 

(1) Ric ft + h = 0{x2) uniformly on the closure of this face, and more 
precisely Ric h + h G a^A"; 

(2) when x\ — > 0, then h — ho G xiA 2 ' a , where ho is the model metric 
\23\) given by the formula 

'dx\ 7l / dx 2 , 72 Jjf_ 

X 2 + 4^2 ^ 4 x f) ±A2 \ x 2 + 4x 2 + 4x 4 

Proof. This is a direct consequence of the formulas (f30|) and (j3lT) for Ric 9 , 
since the solution (91,72) depends smoothly on the parameter q2 G X2. □ 



We conclude this section with a comment about regularity. Exactly as in 
Remark [22] at the very end of § HI the solutions obtained in Theorem [33] are 
poly homogeneous at the boundaries of the codimension one faces provided 
the (Gi x G2)-conformal infinity data on X is smooth. The proof is identical 
to the one for the uncoupled AKH Einstein equations. 

6.2. Strongly asymptotically Einstein metrics. We are now ready to 
define, given any small deformation c of the given G\ x GVconformal struc- 
ture c° on X, a global, approximately Einstein metric g on M. By The- 
orem [33] and Corollary 134] we extend the data (71,72) on the two faces 
F\ = M\ x X2 and F2 = X\ x M2 to get pairs (91,72) and (71,92) solving 
the system (I32p on each face. Therefore, the two metrics on M, 



h \ ^ 2 ^ \ . X f dX i , 71 , Vl \ , 

hl=9l + X2 {^ + 4x 2 + 44)> /l2 = Al vT + ^! + MJ +52 
defined in the neighbourhoods X2 <C 1 and 27 <C 1 of iq and F%, respectively, 
satisfy 



hi - hoe xiA 2 ' a , h 2 - h G x 2 A 2 ' a , 



where 



dx\ 71 *7i V . / <fa| 1 72 T?| 
x 2 + 4x 2 + Ax\) + 2 V x 2 + 4x 2 + 4x 4 2 
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is the initial model metric (|23f) . It remains to glue h± and /12 in the region 
{x\ -C 1,X2 <C 1} where they both exist. Choose some cut-off function \ 
such that x( x ) = 1 f° r x < 1/2 and = f° r x > 2, and consider the 
metric 

* = x£)fa + (i-x(S))h, 

now defined in a neighbourhood of the whole boundary F\ U F2. In the 
region 1/2 < x%/x\ < 2, all the derivatives of x( x 2/xi) remain bounded for 

the metric + hence in that region h — ho £ xiA 2 ' a (or equivalently 

x \ x 2 

a^A 2 '"), and therefore Kic h + h = 0(x\), or more precisely Ric^+ft. G x\K 2,a . 
Globally, in a neighbourhood of F\ UF2, 

(33) Ric' 1 + /» € x£a£A a , i/ = i. 

We now generalize this model. Fix < v < 1/2. 

Definition 35. ^4 metric g on M is strongly asymptotically Einstein if it 

differs from the metric h defined above by a term in x\x^K 2,,a . 

In particular, the Ricci curvature of any such metric satisfies (|33|) . Note 
too that any other reasonably method of patching h\ and h% together near 
the corner yields a metric b! which is strongly asymptotically Einstein in 
this same sense. 

7. Generalized Laplacians on near-product hyperbolic spaces 

We now discuss the construction of a parametrix for L 9 when g is a 
strongly asymptotically Einstein perturbation of a product hyperbolic met- 
ric. Our goal is to show that L 9 is invertible between two weighted Holder 
spaces. We do this in the following steps. First, we analyze the Schwartz 
kernel of the inverse of L 9 when (M, g) is exactly product hyperbolic using 
a contour integral representation; we go on to obtain conormal bounds for 
this Schwartz kernel on the 'product hyperbolic double space' M^ h . This 
serves as an ansatz for the parametrix of L 9 when the metric g is weakly 
asymptotically product hyperbolic. We introduce a stronger condition on g 
of being near product hyperbolic, and show that under this hypothesis we 
can construct a parametrix with Schwartz kernel conormal on Mp h which is 
an inverse of L 9 up to a compact error term. The final step is to show that 
each of these operators are bounded between weighted Holder spaces, which 
implies that L 9 is Fredholm on these spaces. The fact that when g is a prod- 
uct we have an exact inverse for L 9 on L 2 which is bounded between these 
Holder spaces shows that L 9 is invertible between these spaces then. The 
parametrix construction varies continuously with g, so we conclude that L 9 
remains invertible when g is near-product hyperbolic and sufficiently close 
to a product metric. 
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7.1. The inverse of L when g is a product. We begin with an examina- 
tion of the structure of the inverse of L 9 on L 2 (M) when (M, g) is a product 
of AKH spaces. This is mostly a review of the analysis in [21] . We first 
present a contour integral representation for the L 2 inverse G of L 9 involv- 
ing the resolvent families of the operators L 9i on each factor; from this we 
deduce estimates for the pointwise off-diagonal behaviour of the Schwartz 
kernel of G using the analogous estimates for the Schwartz kernels of the 
two constituent resolvents. We do not make an effort to obtain the most 
precise pointwise estimates on G here, but see [21] and [14] for more on this. 

A representation formula for the Green function. Let (M, g) be a 

product of two AKH Einstein spaces, with linearized gauged Einstein oper- 

o 

ator L 9 = V*V — 2R, acting on sections of 

Sym 2 T*M = Sym 2 T*M 1 © (T*Mi ® T*M 2 ) © Sym 2 T*M 2 . 
The operator L 9 preserves the three summands and acts by 

1,91 _|_ (y M 2y^M 2 /yAfi^*yMi _|_ ^jM 2 y^jM 2 ^yMi\*yMi _|_ jjj 2 

respectively. In each of these three cases, it has the form 

L 9 = Li I B2 + I Bl ® L 2 , 

where Lj acts on a Banach space Bi of sections of a bundle Ei on Mi, and 
L 9 acts on the completed tensor product B\ ®B 2 of sections of E = E% ® E 2 
on M. 

For the moment, let Bi = I? {Mi, Ef, dV gi ). Denote by Ri(fi) the resolvent 
family (Li — This is a holomorphic family of bounded operators on Bi 

for [i in the resolvent set C \ spec (Li); according to Theorem 1141 

spec (Li) = {Ay}^i U [oti, oo) 

for some aj > 0, with A^- G K, Ajj < a^. Nondegeneracy of (Mi,gi) for the 
Einstein problem is the assumption that ^ {Ajj}. 

The resolvent family of the product operator, R(\) = (L — A) -1 , can be 
expressed as a sort of convolution of the resolvents on the two factors. More 
precisely, it is proved in [21] that 

(34) R( A) = - -L / Ri(fi) i?2(A - n) dfi, 

where T\ is a contour lying in the common region of holomorphy of the 
two factors in the integrand with ends converging linearly to ±ioo and such 
that the spectrum of L\ lies entirely on one side and the spectrum of L 2 lies 
entirely on the other side. For us it suffices to take A = 0. In the simplest 
situation, neither L\ nor L 2 have any negative eigenvalues, and in this case 
we take Tq = iM.. The general case, where one or the other does have such 
eigenvalues, requires a slight modification to this formula. 

Setting R(0) = L~ l , and G its Schwartz kernel, if / G ^§°(M), then 
u = Gf is the unique L 2 solution to the equation Lu = f. 
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Estimates on G. Fix local coordinates Zj = (xj, yj) near a boundary point 
of Mj, and denote by Kj(zj, Zp fi) the Schwartz kernel of Rj(/J-). Also, re- 
place fi by ifx. Then (assuming neither operator Lj has negative eigenvalues), 
the Schwartz kernel of G equals 

1 f°° i 

(35) G(z 1 ,z 2 ,z' 1 ,z 2 ) = — / K 1 (z 1 ,z' 1 ,ifi) K 2 (z2,z' 2 ,-ifj,)dfi. 

We now obtain pointwise estimates for G using the bounds on Kj(zj, z'j,ifj,) 
in Proposition [HI Our goal is to prove that G lifts to a conormal distribution 
on a certain blowup of M 2 which is conormal at all boundaries, polyhomo- 
geneous along the lifted diagonal, and has a leading polyhomogeneous term 
at the front faces with conormal remainder. 

First observe that the integral (|35p converges in the Banach space of 
bounded operators on L 2 (M;dV g ). This follows from the elementary esti- 
mate 

ll-^i(^)IU(i 2 (Mj)) - r-f) 
\t L \ 

which is a direct consequence of the spectral theorem for the selfadjoint 
operator Lj. This already uniquely specifies G as an element of &(M 2 ). 
Next, introduce the product-hyperbolic double space 

M p \ = (Mi)| ie x (M 2 )i 20 ; 

this manifold with corners is simply the product of the Kj-0 double spaces 
of the two factors. It has six codimension one boundary faces: the two front 
faces, which are the boundary hypersurfaces intersecting the lifted diagonal 
diag ph : 

ffl = ff((Mi)| l6 ) x (M 2 )| 20 and ff 2 = (Mi)| l9 X ff((M 2 )| 26 ), 

and the four side faces, which are the products of the side faces of one factor 
with the interior of the other factor. We denote the side faces by -Bioj and 
Boiji j = 1)2, and for uniformity also write ffj = Buj; the j signifies 
that the face in question comes from a boundary face in the Mj factor. 
Defining functions for any one of these faces will be written p pq j, j = 1,2 
and pq = 10,01,11. 

We shall prove that G is conormal at all boundary faces of M 2 h and along 
the lifted diagonal. The conormal estimates here are not quite sharp; there 
is a subtle cancellation in (j35[) . explored more carefully in [21], which leads 
to vanishing rates which are slightly better (by a logarithmic factor), but 
this is not needed here. 

Proposition 36. The lift of the Schwartz kernel G to the space M 2 h is 
an element of £/^(M 2 h , diag ph ), where a = (c pg ,j) is the index set with 
a w,j = a oi,j = n j + dj — 1, ciij = 0, and where ff = ffi U ff 2 . 

Proof. To separate out the contributions from the near-diagonal parts of 
each factor, we use standard results concerning the symbol calculus with 
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spectral parameter, cf. [26]. Write Kj = Aj + Bj where Aj contains the full 
diagonal singularity, is supported near the lifted diagonal in (-Mj)jg e , and is 

smooth across the front face of this double-space, while Bj £ ((Mj)^ e ), 

a . = {{nj + d . _ i )/2 + ( nj + dj - l)/2 + ^'). 

The integral of ^4i(i/i)^42(— i//) is a distribution supported in a neighbour- 
hood of diag ph which intersects only the front faces of M ph but not the side 
faces. It may be estimated directly (just as for the analogous computation on 
a product of two compact manifolds) using the oscillatory integral represen- 
tations in the conormal bundle of diag ph , and the result is clearly smooth 
up to the front faces. This is the Schwartz kernel of a pseudodifferential 
operator of order —2, as required. 

Next, the norm of Aj as a map on any fixed weighted Sobolev or weighted 
Holder space decays like l/|/i|, cf. [26]. (Note that this is equivalent to the 
situation in the compact case because of the support properties of Aj.) From 
this we deduce that the integrals of Ai(ifi)B2(—ifi) and A2(in)Bi(—i/j 1 ) are 
smooth on the interior of M ph . By taking any number of derivatives with 
respect to b- vector fields on the B factor, we obtain the correct conormal 
estimates too. 

As for the integral of Bi(ifi)B2(— i/i), divide the contour Tq into a compact 
portion, where < /io, and its remaining noncompact ends > /j,q. The 
conormal estimates on Bj hold locally uniformly in /i, so in particular 

/ B 1 (zi,z' 1 ,ifj,)B 2 {z 2 ,z 2 ,-ifj,)dfj, £ ^ n (M ph ). 
J\^\<lio 

On the other hand, using the exponential bounds on Kj from (I16p we imme- 
diately deduce that the integral over |^| > fj,Q satisfies the same ^° bound. 
The estimates for higher tangential derivatives required to check conormal- 
ity are obtained in exactly the same way, using the corresponding pointwise 
bounds for the higher tangential derivatives on each of the two factors. 
This completes the estimation of the Schwartz kernel G. □ 



Modifications when either factor has negative point spectrum. If 

infspec(Li) or infspec(Z/2) are negative, then (|35h needs to be altered 
slightly. To understand this, begin by noting that (|34p remains valid when 
A = ie and the contour separates the sets ie — spec (Li) and spec(L2) m two 
different half-planes. In fact, we take the contour to be the union of the 
vertical rays [2ie, ioo), (— ioo, — ie] and two long thin half-ellipses, one in the 
first quadrant with minor axis connecting ^ie and 2ie, and the other in the 
third quadrant, with minor axis connecting —\it and ^ie. Now let e \ 0; 
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we arrive at the formula that (L\ + L?) 1 is equal to the sum 

- — / R 1 (-n)R 2 {^)dn 
2 ™ J_ ioo 

Here //^ are the (negative) eigenvalues of Lj, and Pj(/z^) is the orthogonal 
projection onto the corresponding eigenspace in the i th factor. 

The main result in the last subsection, that the Schwartz kernel of (L 9 )^ 1 
is conormal on M 2 h , clearly remains valid. The analysis of the main term, 
which is the first summand in this expression, is exactly the same as before. 
The remaining terms are much simpler to analyze since each is a simple 
tensor product of one term which is polyhomogeneous on one factor and a 
finite rank polyhomogeneous term on the other factor. 



7.2. The parametrix in the near-product case. If g = g\ + g<i is a non- 
degenerate product hyperbolic metric, we have established that the Schwartz 
kernel of the L 2 inverse of L 9 has only conormal singularities on the prod- 
uct hyperbolic double space. We now show that a similar structure theorem 
holds for a parametrix for (L 9 )^ 1 , i.e. an approximate inverse up to compact 
error terms, when g satisfies a condition we shall call being 'near product- 
hyperbolic'; the precise definition is given later. This parametrix H will 
have Schwartz kernel in ^ p ^(M) + =*#ff(M^ h ). The construction of H uses 

the structure of (L 9 )^ 1 in the product case as an ansatz and a tool. We 
describe this now, deferring the statement of the final result to a theorem 
at the end of this subsection. 

The procedure is very much the same as in the AKH case. The Schwartz 
kernels of general ph pseudodifferential operators are defined to be distribu- 
tions on Mp h which are conormal at the boundaries and with a polyhomoge- 
neous expansion at the diagonal. There is a small calculus of operators with 
Schwartz kernels supported near the lifted diagonal, and a large calculus 
which also admits operators with Schwartz kernels conormal up to all faces, 
and with positive vanishing order at the front faces. The corresponding 
decomposition of H is written H\ + H 2 . 

Let g be a weakly asymptotically product hyperbolic metric on M. The 
initial approximation H\ to (L 9 )" 1 is obtained via the standard elliptic 
parametrix construction using the symbol calculus on the conormal bundle 
of diag ph in M 2 h . This uses the uniform invertibility of the ph-symbol of 
L 9 up to the two front faces, which is valid for any weakly asymptotically 
product hyperbolic metric. If g is only polyhomogeneous at the boundaries 
of M, then Hi G ^ (M^ h , diag ph ), where t W j = r i,j = oo and mj = 0, 
J = 1,2. 
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Set L 9 H\ = I — Q. By construction, Q E sd^ (M? h ). The correction term 
i?2 is chosen so that L(-ffi + H2) = I — Qi, or equivalently, LH2 = Q — Qi, 
where H 2 E (M^, diag ph ) and Q x E ^(M p \); here a UJ = 0, or UJ > 0, 
and all other 0p 9 j = <Tpg,j = + — 1 for pq 7^ 11. The vanishing of the 
Schwartz kernel of Q\ to some positive order at the front faces yields its 
compactness on weighted Holder spaces. 

To determine H2, restrict the equation L 9 LH 2 = Q ~ Qi to each of the 
two front faces. Set 

N^W) = L%. , NjiQo) = Q | ffj , 3 = 1, 2. 

Then we must solve the two equations Nj(L)Hj = Nj(Qo). The right hand 
side is a smooth function on the interior of Sj; it vanishes to infinite order 
at the intersection with the side faces, and is smooth up to the other front 
face, or at least, the dependence in this direction is exactly as regular as the 
metric g near the corner. 

Let us begin by analyzing the structure of Nj(L). Fixing j = 2 to be 
definite, the front face fT 2 is a product (-Mi)jg- e x ff (M2). The second factor, 
the front face of (M 2 )^ q, is a fibration with base space the diagonal of 
(OM2) 2 and each fibre naturally identified with the hyperbolic space K2H m2 . 
The lift L 9 from the left factor of M to M p \ acts on the 'left factor' of Mj in 
each KjQ double space. Its restriction to ff2, which makes sense since it acts 
tangentially to that boundary, is a sum of derivatives of two types: some act 
on the left factor of M\ in (Mi)^ i0 and the others act on the K^-ff" 12 fibres 
of ff(M2). In particular, the dependence on all other variables - namely, in 
the right factor of M\ and in the diagonal of (OM2) 2 - is purely parametric. 
The key assumption is that N 2 (L 9 ) = £2,1 + L 2} 2, where £2,1 is an operator 
acting on the left factor of M\ (with coefficients depending parametrically 
on the right factor of M\ and the diagonal of (SM2) 2 ) and -£^2 ^ s the gauged 
linearized Einstein operator or the rough Laplacian on W>2H m2 . Similarly, 
N\{L 9 ) = L± t i + Li t 2, where these summands have analogous properties. 

Let us investigate when this condition holds. Using interior coordinates 
Z\ on M\ and coordinates z 2 = (#2,2/2) near the boundary on M 2 , we have 
near #2 = 0, 

L 9 = 5>W(* ljlte )5£ +Y J ^ n (z 1 ,y 2 )(x2d X2 y(x2Yf(xlzr' + ff(x 2 ); 

where (Yi, . . . , Yd 2 (m 2 -i)) an d (Zi, ■ • • , ^d 2 -i) are Dases of sections of & 2 a nd 
its complement, respectively, in TX 2 , and we are using obvious multi-index 
notation. Notice that we are using the product structure Mi x dM 2 for 
this boundary face. Near product hyperbolicity is simply the condition that 
there are no cross-terms at X2 = 0, and similarly at x\ = 0. Note that by 
naturality of generalized Laplacians, the second summand here, at #2 = 0, 
must be the linearized gauged Einstein operator on K^-ff" 12 - 

Lemma 37. Let g be a strongly asymptotically product hyperbolic metric, 
as constructed in § \ 6.SX Then g is near product hyperbolic. 
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Proof. By the global rigidity assumption for the conformal infinity data, we 
are fixing the product structure on each boundary hypersurface. Consider 
the construction of g at Mi x dM<i- The preliminary extensions of the 
conformal infinity data (71,771) and 72/^2 + r ?2/^2' regarded as metrics on the 
bundle TM\ (B K2@ Tqm 2 M2 over this face, make these subbundles orthogonal 
to one another. The correction terms h and <p are chosen to preserve this 
orthogonality. □ 

Now, suppose that L g has this near product form. Consider its lift to M^ h ; 

near ff 2 we replace (22,22) with projective coordinates, so 2/2 = 2/2 + G(x-i) 

and the derivatives with respect to X2, Y% and Zj in the second summand all 

become tangential to K^-ff" 12 . Setting X2 = gives N2(L). The coefficients 

da then depend on [z\, 2/2), but not the second hyperbolic space factor; the 
(2) 

coefficients a - L potentially depend on all variables on this face, although 
the derivatives just act tangentially to the hyperbolic spaces, but 1/2,2 is 
simply the linearized gauged Einstein operator or the rough Laplacian on 
!€2-ff m2 , and hence is independent of z\ and 2/2- in other words, if g is any 
near product hyperbolic metric, in particular a strongly asymptotically Ein- 
stein metric from Theorem [331 then the restriction of Nj(L 9 ) to the fibres of 
the front face ff j is of product type with respect to the product structure on 
each fibre. A corresponding statement is true for any generalized Laplacian 
P 9 associated to g. 

Having established this, we now proceed to solve the two normal equations 
using the results of § 17.11 More specifically, since Nj(L g ) is of product type, 
the results of § 17.11 show that Nj(H2) = N j(L 9 ) N j(Qo) is conormal on ff,- 
and smooth up to the corner ffiHfi^. To verify this last statement, note that 
the product analysis of § 17.11 gives conormality on each slice M\ x M.2H m2 
(when j = 2, for example), but the dependence on the other variables, i.e. 
the 'right' factor of Mi in (-Mi)jjL e and the base space variable in ff((M 2 )jg @, 
are parametric, and the solution is as smooth as the dependence on these 
parameters. As the variable in the right factor of Mi converges to the 
boundary, the space on which N2{L 9 ) acts converges locally uniformly to the 
product of the two hyperbolic spaces; the uniform estimates on iV^Qo) and 
on N2(L 9 )~ 1 show that the solutions converge smoothly too, and that the 
limiting solution is conormal on the product of the radial compactifications 
of the hyperbolic spaces. In particular, the corner ffi n ff 2 is a bundle over 
dM\ x DM2 with fibre the product of hyperbolic spaces, and the limit of the 
solutions from either side, ff 1 or ff 2 are the same. Therefore 

^l(^2)| ffinff2 = iV 2 (if 2 )| ffinff2 . 

This argument proves the existence of conormal functions on the two front 
faces which are compatible at the corner. Hence we may extend these to a 
Schwartz kernel H2 in the interior of M which is conormal to all boundary 
faces, H2 G ^(M^), and which restricts to iV^i^) at ffj. 
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We have now constructed H = H\ + H 2 which satisfies L 9 H = I — Qi, 
Qi G £/ a (Mp h ). A left parametrix is obtained by taking adjoints. This 
completes the parametrix construction when g is near product hyperbolic 
and polyhomogeneous. 

As explained earlier, we can certainly restrict to studying polyhomoge- 
neous metrics with smooth conformal infinity data but we comment briefly 
on how to extend this proof to the case where g is the sum of a polyho- 
mogeneous AKH metric and a perturbation term k G x^x^A^q, where 
< Uj < rij + dj — 1. (Note that this is the regularity for the strongly 
asymptotically Einstein metrics.) 

The symbol calculus step goes through immediately for symbols with this 
regularity. Indeed, the perturbation term k appears only in the parametric 
dependence along the diagonal and does not occur in the leading terms near 
the front faces, so Hi decomposes as a principal polyhomogeneous term and 
another, vanishing to some positive order at these faces but which is only 
bounded in Holder norm. The restriction to the front faces of the error 
term Q = I — LH\ does not depend on these lower order terms, and the 
normal operators Nj(L) are independent of them too, except as above in 
their dependence on parameters, which means that the second step carries 
through exactly as before. 

We have proved the 

Proposition 38. Let g be a strongly asymptotically Einstein metric as con- 
structed in §EH Then there is a parametrix H = Hi + H 2 for L 9 , with Hi G 
sf£ (M^ h ,diag ph ) and H 2 G -^s( M ph)- The error term Qi = I-L(H 1 + H 2 ) 
isin^i(Ml h ). 

7.3. Function spaces and mapping properties. We now show that the 
pseudodifferential inverses or parametrices for L 9 , when g is product, or near 
product, hyperbolic, are bounded on weighted 'geometric' Holder spaces 
associated to the metric g. 

The definition of these Holder spaces is essentially identical to the one in 
§ O the only difference when the weight parameters are being that to define 
A^ when k > 0, we allow differentiations by arbitrary smooth vector fields 
which are locally given as finite combinations of smooth products of elements 
of '^k 1 e(Mi) x ^K 2 e(M2). If v\ and v 2 are any two weight parameters, we 
also define x" 1 x^ 2 (M) in the obvious way. 

If Ki = H, then the Holder spaces corresponding to two different near 
product hyperbolic metrics g and g' are not quasi-isometric to one another, 
unless the corresponding distributions Q>\ and S>[ on dM± x M 2 are diffeo- 
morphically equivalent. 

Proposition 39. Let H be an operator which has its Schwartz kernel in 
=G^f (Mp h , diag pll ), where a is the usual weight family with anj = and all 

other ap q j = (rij + dj — l)/2 + 8 J , and with pseudodifferential order —£ 
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along diag ph . Let < i/j < (jij + dj — l)/2 + 5 3 . Then 

H : x?x?A k p £(M) — > x?x?A k +*> a (M) 

defines a bounded map. If Q has Schwartz kernel in £f a (M 2 il ) where &ii t j > 
but all other weights are the same as for a, then Q defines a compact 
mapping between these same spaces. 

Proof. As before, write H = Hi + Hi- The boundedness of Hi follows 
from the standard boundedness of pseudodifferential operators of on Holder 
spaces as well as the fact that the support conditions on this Schwartz kernel 
means that this part of H only spreads supports by a fixed amount. 

The proof of the < ^ , ° bound for H%f when / is in this weighted Holder 
space proceeds by a direct and elementary estimation of this integral. Any 
higher derivative of H^u with respect to a Ki0 vector field on the first factor 
or a vector field on the second factor, or an iterated combination of 
such vector fields, is handled by simply noting that any of these vector fields 
applied to H2 gives a Schwartz kernel of exactly the same form, with the 
same orders of vanishing at all the side faces. 

Finally, to prove that Q is a compact operator, the extra vanishing at the 
front faces gives that 

Q : x^Aj^M) — > x^x^A^iM) 

for any si,S2 £ K. This range space on the right includes compactly into 
the domain space on the left (provided S2 > si), using the Arzela-Ascoli 
theorem. □ 

Theorem 40. Let g be a near product hyperbolic metric which is sufficiently 
close to a nondegenerate product hyperbolic metric go on Mi x M2. Then 

L 9 : x^x^A 2 ^{M,S 2 T*M) — ► xf^ 2 A°^(M, S 2 T*M) 

is an isomorphism. 

Proof. If g is any near product hyperbolic metric, then we have constructed 
a parametrix H 9 for L 9 so that L 9 H 9 = H 9 L 9 = I — Q 9 . We have also 
proved that H 9 is a bounded operator between the space on the right and 
the space on the left above, and that Q 9 is compact between these same 
spaces. This proves that the map L 9 is Fredholm. 

The construction of H 9 depends continuously on the metric (with respect 
to some sufficiently strong topology). Since Q 9 = when g = go, we can 
make the norm of Q 9 as small as desired when g is sufficiently close to go, 
which implies that L 9 is invertible. □ 

8. Solving for the Einstein metric 

The remainder of the proof of the Main Theorem proceeds very much as 
in the analogous arguments in § |4] and § [HJ 
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Let (M = Mi x M2, g = gi + 52) be a product of AKH Einstein metrics. 
Let L 9 and L 9i be the linearized gauged Einstein operators for g and the 
two component metrics g^. We assume that 

i spec (L 9 ) U spec (L 91 ) U spec(L 92 ), 

and that at least one of K4 or K2 is not equal to EI or O. Denote by c the 
{G\ x G2)-conformal infinity data on X = X\ x X 2 = dM\ x dM 2 . 

Theorem 41. Under all these conditions, let c' be any globally integrable 
(G± x G2)-conjormal infinity data on X which is ff 00 and sufficiently close 
(in c to 2,a norm) to c. Then there is a unique near product hyperbolic Ein- 
stein metric g' which is close to the near product hyperbolic asymptotically 
Einstein metric g' with conformal infinity data c' constructed in § \ 6.2\ . 

Proof. Let g' be the strongly asymptotically Einstein metric constructed in 
§ 16.21 with conformal infinity data c'. Write g' = g' + k, where 

kex^ 2 A 2 ^(M,Sjm 2 (T*M)). 

Write the gauged Einstein equation as 

N' 9 '(k) = mc(g' + k) + (g' + k) + {5~ 9 ' +k )*B 9 \k) = 0. 

The linearization at k = is the generalized Laplacian L 9 ' . 
We have that 

N 9 ' : x?x?A%£(M, Sym 2 (T*M)) — » x u 2 2 A° p £ (M , Sym 2 (T*M)) 

is a ff 1 map from a neighbourhood of in the domain space, and that L 9 ' 
is a bounded linear map between these same two spaces. 

According to Theorem 1401 if c' is sufficiently close to c, this linearization 
is an isomorphism. Furthermore, the norm of its inverse is bounded away 
from zero, uniformly as c' — > c. The inverse function theorem implies that 
there is a unique solution k to N 9 (k) = with k near 0. □ 

Appendix A. Osculating quaternionic coordinates 

In this brief appendix we prove Lemma [6j This uses an idea close to that 
used to find normal coordinates in Riemannian geometry, and should be a 
general fact for all the so called "parabolic geometries". 

Fix a quaternionic contact structure on y 4m_1 and any metric in this con- 
formal class on the distribution £F. We use the Tanaka- Webster type con- 
nection V from [21 chapter 2] , see also [U [7] for the special case of dimension 
7. The contact distribution 3) is the kernel of three 1-forms (7/1, 772, 773), and 
has a privileged supplementary subspace generated by three "Reeb vectors 
fields" R\, i?2, -R3 which are uniquely specified by the conditions r]i(Rj) = 5ij 
and (ij^drjj + ir drji)\@ = (in dimension 7 a weaker condition is placed on 
these). This connection satisfies: 

• V preserves the distribution £F and the quaternionic structure on 3>\ 
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• the torsion T v of two horizontal vectors X, Y £ <2l is given by 

3 

T% x = J2 d Vk(X,Y)R k . 
1 

However, since these conditions only place restrictions on the derivatives in 
horizontal directions, the connection is not unique and different extensions 
are possible. 

Now fix a point p £ Y. We may assume that the Ri, and hence also the 
rji, are parallel at p, i.e. 

VR l (p) = 0, i = 1,2,3. 
The connection determines the exponential map in horizontal directions, 

exp p : % — >Y, 

by solving the differential equation VcC = 0, c(0) = p, c(0) = X, and setting 
exp p (X) = c(l). Linear coordinates (y\, . . . , y^m-i) on give a coordinate 
system on the image S of a small ball by exp . Then, at p one has 

and hence 

3 

l 

In particular, still at p, 

d yiVk(d yj ) = m(^dy t d yj ) = ^dr] k (d yi ,d yj ), 
from which we deduce that for y £ S, the projection Y{ of d Vi on 3) satisfies 

1 3 

Y 3 = d v 3 - jE^fti.^i)^ + °(\y\ 2 )- 

1 

This can be interpreted as saying that if we write the standard quaternionic 
contact structure Go in coordinates (aj,yi) on the Heisenberg group, and 
denote by Y® the standard horizontal vector fields as in Q and (fTUj) . then 
along S one has 

rj = e + O(\y\ 2 ), Y i = Y? + 0{\y\ 2 ). 

Choosing transverse coordinates (a\, cr 2 , 03) so that d ai = Ri along S, we 
get the same result in a neighborhood of p with an error term 0(|y| 2 + |cr|). 
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